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STOCHASTIC  CROSS-SECTIONS  BASED  ON  THE 
SMALL  SLOPE  APPROXIMATION:  THEORY 


1  Introduction 

1. 1  Background  and  motivation 

The  reverberation  of  underwater  sound  is  frequently  modeled  using  ray  theory. 
Particularly  in  shallow  water,  such  a  scenario  involves  bounces  from  stochastic  rough 
surfaces  such  as  the  air-sea  interface  and  the  ocean  bottom.  The  stochastic  incoherent 
scattering  cross-section  as  a  function  of  the  incoming  grazing  angle  and  the  scattered 
angle  is  used  to  specify  the  intensity  of  the  field  that  is  incoherently  scattered  by  the 
surface  roughness. 

For  decades,  perturbation  theory  based  on  small  surface  height  approximation  [1]  has 
been  used  to  calculate  this  stochastic  incoherent  cross-section.  This  result  depends  on 
surface  roughness  on  the  order  of  the  Bragg  scale.  The  Bragg  scale  corresponds  to  the 
spacing  of  a  diffraction  gating  that  produces  primary  maxima  in  the  direction  specified. 
Such  scattering  is  known  as  Bragg  scattering.  In  general,  the  effective  diffraction 
gratings  are  tilted  by  larger-scale  features,  but  perturbation  theory  ignores  this.  However, 
at  high  grazing  angles,  this  effect  can  be  significant  after  only  a  single  scattering  event  [2]. 
Consequently,  there  is  often  a  need  to  improve  the  accuracy  of  the  cross-section  by 
incorporating  the  effects  associated  with  the  large-scale  tilts  that  modulate  Bragg 
scattering.  The  challenge  is  therefore  to  develop  formulations  of  the  scattering  problem 
that  both  improve  on  perturbation  theory  by  incorporating  large-scale  tilts  and  also  lead 
to  tractable  calculations. 

This  paper  imposes  the  tractability  requirement  set  in  reference  [3]:  no  more  than  a  2- 
dimensional  integration  for  2-dimensional  surfaces.  In  most  cases  considered  below,  this 
requirement  is  exceeded,  because  spherically  symmetric  geometry  is  assumed  and  then 
used  to  reduce  the  two  dimensional  integral  to  a  one-dimensional  integral.  Furthermore, 
spectra  are  chosen  so  that  the  associated  autocorrelation  functions,  which  appear  in  the 
intergrand  of  this  primary  (one-dimensional)  integral,  have  closed-form  solutions  or  at 
worst  involve  one-dimensional  integrals  of  finite  domain.  As  in  reference  [3],  we  will 
emphasize  one  particular  approximation  that  leads  to  expressions  for  the  scattering  cross- 
section  that  are  both  consistent  with  our  tractability  requirement  and  incorporate  large- 
scale  tilt.  This  approximation  is  the  lowest-order  small  slope  approximation.  It  was 
developed  by  Voronovich  in  the  early  1980s  [4,  5].  Since  then,  it  has  become  a  widely 
used  and  well-developed  technique  for  calculating  stochastic  cross-sections  for  scattering 
from  quasi-planar  rough  interfaces  of  the  sort  that  occur,  for  example,  on  the  ocean’s  air- 
sea  interface  [6,  7,  8,  9,  10,  11,  12]  or  ocean  bottom  [13,  14,  15,  16].  (See  Section  I  of 
reference  [6]  for  a  discussion  of  how  the  small  slope  approximation  incorporates  large- 
scale  tilt.) 
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The  small  slope  approximation  for  the  cross-section  consists  of  the  sum  of  a  tenn 
corresponding  to  specular  reflection  and  an  incoherent  scattering  cross-section.  The 
former  is  proportional  to  a  S  -function,  and  the  coefficient  of  this  8  -function 
corresponds  to  the  fraction  of  the  incoming  field’s  intensity  that  is  specularly  (i.e., 
coherently)  reflected  from  the  interface.  This  coefficient  accounts  for  the  energy  that  is 
transferred  from  the  coherent  to  the  incoherent  field,  and  it  must  be  a  fraction  between 
zero  and  one.  This  coefficient  is  relatively  easy  to  evaluate. 

The  incoherent  scattering  cross-section  is  a  function  of  the  incident  angle  and  the  angular 
direction  of  the  observer.  It  is  proportional  to  the  intensity  of  the  field  incoherently 
scattered  into  this  direction,  and  it  takes  the  form  of  a  product  between  a  prefactor  and  an 
integral.  All  infonnation  about  the  nature  of  the  media  involved  and  consequently  the 
boundary  conditions  resides  in  the  prefactor,  which  is  a  straightforward  algebraic 
expression  and  is  consequently  easy  to  evaluate.  The  integral  is  the  same  for  all  media 
and  boundary  conditions,  but  it  contains  information  about  the  surface  roughness. 
Specifically,  for  a  quasi-planar  rough  surface  defined  by  a  surface  height  function  h  ( x) , 
it  is  a  two-dimensional  integral  over  the  horizontal  plane  x  with  the  autocorrelation 
function  ^(x)/?(0)^  appearing  in  the  integrand.  (As  discussed  in  reference  [6], 
translational  invariance  and  Gaussian  statistics  are  implicit  in  this  result.) 

If  the  autocorrelation  function  is  known  in  closed  form,  then  we  have  met  our  tractability 
condition.  Unfortunately,  it  is  rarely  convenient  to  directly  characterize  a  random  rough 
surface  by  explicitly  specifying  the  autocorrelation  function.  The  surface  roughness  is 
more  typically  characterized  using  the  Fourier  transfonn  of  the  autocorrelation  function. 
This  quantity  is  called  the  surface  roughness  spectrum  S(  l_) .  If  the  spectrum  can  be 

transfonned  back  to  give  a  closed  fonn  of  the  autocorrelation  function,  then  our  tractable 
condition  is  automatically  satisfied.  A  central  focus  of  this  paper  will  be  to  develop  and 
examine  spectra  that  are  both  physically  realistic  and  either  meet  or  at  least  nearly  meet 
this  condition. 


1.2  The  tradeoffs  involved 

As  noted  above,  throughout  most  of  this  paper  (Sections  3  through  6),  a  considerable 
simplification  is  obtained  by  invoking  spherical  symmetry  to  reduce  the  two-dimensional 
integral  over  x  to  a  one  dimensional  integral  over  the  magnitude  |x| .  This  integral 

contains  a  very  slowly  decreasing  oscillating  integrand,  and  it  is  still  fairly  difficult  to 
evaluate.  Spectra  must  be  formulated  in  such  a  way  that  the  integrand  is  easy  to  evaluate. 
The  various  spectra  discussed  below  involve  various  tradeoffs,  but  two  are  most 
significant.  The  first  tradeoff  involves  the  ease  of  obtaining  the  autocorrelation  function 
(i.e.,  tractability)  versus  the  ability  of  the  spectrum  to  incorporate  as  many  known 
attributes  of  actual  the  physical  rough  surface  spectrum  as  possible.  Models  whose 
characterization  of  the  surface  is  comparatively  less  faithful  may  still  give  accurate  cross- 
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sections  under  many  important  scenarios,  but  they  are  unlikely  to  be  as  broadly 
applicable  as  those  that  incorporate  a  more  accurate  and/or  detailed  description  of  the 
interface.  Thus,  the  most  significant  tradeoff  between  spectra  involves  tractability  versus 
applicability  to  a  broad  variety  of  physical  scenarios.  The  second  tradeoff  will  involve 
spectra  that  offer  comparable  tractability,  but  employ  different  simplifying  assumptions. 
In  this  case,  the  goal  will  be  to  match  a  given  application  with  the  optimal  model.  Even 
in  this  context,  a  premium  will  be  placed  on  models  that  have  the  broadest  applicability. 

An  example  of  the  first  tradeoff  was  explored  in  reference  [6],  A  pair  of  spectra  were 
considered  and  the  techniques  for  solving  the  associated  integrals  in  the  cross-section 
were  developed.  This  reference  constructed  two  spectra  that  satisfied  certain 
characteristics  assumed  to  hold  for  the  surface-roughness  spectrum  of  the  air-sea 
interface.  In  the  most  tractable  scenario,  it  was  assumed  that  only  the  (transversely) 
isotropic  tail  of  the  spectrum  was  relevant.  The  integration  was  solved  and  the  associated 
cross-section  involved  no  numerical  integrations.  In  a  second  scenario,  information 
about  the  spectral  peak  as  well  as  the  tail  of  the  spectrum  was  used  to  generate  an 
isotropic  “difference  spectrum”  formed  by  taking  the  difference  between  two  power  laws. 
In  this  case,  the  integrand  in  the  cross-section’s  one-dimensional  integral  from  zero  to 
infinity  appeared  in  closed  form  and  the  intergral  was  solved  using  specialized  methods. 
This  paper  applies  the  pure  power-law  and  difference  spectra  to  parameter  sets  other  than 
those  considered  in  reference  [6],  and  it  furthennore  introduces  classes  of  spectra  not 
considered  in  reference  [6]  (see  Fig.  1). 

While  the  difference  spectrum  can  accommodate  an  arbitrary  peak  location  and  power- 
law  tail  behavior,  there  are  limits  to  the  sharpness  of  the  peak.  The  difference  spectrum 
peak  is  therefore  sharpened  either  by  replacing  it  at  low  wavenumbers  with  a  decaying 
exponential  function,  or  by  replacing  it  entirely  with  a  Pierson-Moskowitz  spectrum 

modified  so  that  the  exponential  e ~a^  is  replaced  with  e  u// . 

In  the  fonner  case,  the  low-wavenumber  exponential  function  generates  an 
autocorrelation  function  with  a  relatively  easy-to-evaluate  one-dimensional  numerical 
integral  (that  is  in  turn  nested  in  the  integrand  of  the  infinite  integral  that  appears  in  the 
cross-section).  Thus,  for  the  cutoff  difference  spectra,  we  have,  relative  to  the  difference 
spectrum,  a  modest  decrease  in  tractability  and  a  modest  increase  in  the  completeness  of 
the  physical  description. 

This  modification  to  Pierson-Moskowitz  spectrum  is  made  so  that  the  autocorrelation 
function  has  a  closed-form  solution.  This  spectrum  therefore  turns  out  to  be  as  tractable 
as  the  difference  spectrum  and  it  also  incorporates  sharper  spectral  peaks  than  any  of  the 
other  spectra  considered  (see  Fig.  1),  but  it  only  meets  these  tractability  conditions  for 
integer  power  laws.  Furthennore,  the  height  of  the  spectral  peak  and  its  location  are 
correlated,  and  so  the  modified  Pierson-Moskowitz  spectrum  is  the  least  flexible  of  all  the 
spectra  considered.  The  primary  tradeoff  between  tractability  and  broadness  of 
applicability  cannot  be  used  to  distinguish  between  the  difference  spectrum  and  the 
modified  Pierson-Moskowitz  spectrum.  These  spectra  accurately  incorporate  different 
aspects  of  reality.  The  difference  spectrum  accurately  models  the  behavior  of  a  wide 
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variety  of  spectral  tails,  while  it  makes  approximations  concerning  the  nature  of  the  peak. 
The  modified  Pierson-Moskowitz  spectrum  will  accurately  reflect  the  shape  of  the  peak, 
while  it  will  be  forced  to  make  approximations  about  the  tail  (say  be  restricted  to  integer 
power  laws).  A  central  theme  of  this  paper  will  be  to  evaluate  this  tradeoff,  and  a  key 
conclusion  will  be  that  the  difference  spectrum  that  accurately  reproduces  a  variety  of  tail 
behaviors  provides  a  more  useful  compromise  than  the  modified  Pierson-Moskowitz 
spectrum  that  emphasizes  peak  shape  at  the  expense  of  describing  only  integer  power-law 
tails. 


Fig.  1  -  The  various  types  of  isotropic  spectra  considered  in  this  paper.  Note  that  this  graph  uses  linear 
scales  for  both  axes.  The  pure  power  law  spectrum  is  given  by  the  dashed  line.  The  difference  spectrum  is 
given  by  the  solid  line.  The  difference  spectrum  is  made  narrower  in  one  of  two  ways;  by  introducing 
either  a  step  cutoff  at  wavenumber  A  ,  or  an  exponential  cutoff.  The  autocorrelation  function 
corresponding  to  the  step  cutoff  is  proportional  to  the  autocorrelation  function  corresponding  to  the 
difference  function  minus  the  contribution  from  the  cross-hatched  and  black  areas.  The  autocorrelation 
function  corresponding  to  the  exponential  cutoff  spectrum  differs  from  that  for  the  step  cutoff  in  that  the 
area  in  black  (i.e.,  under  the  exponential)  contributes  to  the  Fourier  transform  of  the  spectrum.  In  this 
illustration  C  <  0  ,  and  the  exponential  cutoff  spectrum  is  concave  from  the  left.  In  order  to  produce  a  peak 
that  is  both  narrower  and  taller  (relative  to  the  power-law  tail)  than  the  difference  spectrum,  a  modified 
Pierson-Moskowitz  spectrum  is  considered  (dash-dot  line).  This  spectrum  only  meets  our  tractability 
conditions  for  integer  power  laws. 


In  the  final  scenario  considered  here,  the  difference  spectrum  is  generalized  to  include  an 
azimuthal  dependence,  and  the  associated  scattering  cross-section  is  shown  to  be 
reducible  to  a  sum  of  tractable  one-dimensional  integrals. 

The  newly  expanded  technique  has  recently  been  used  to  develop  improved  algorithms 
for  scattering  from  the  air-sea  interface  [17,  18],  and  to  study  the  scattering  behavior  of 
the  ocean  bottom  [2,  19,  20],  An  overview  of  some  of  the  material  discussed  below  is 
published  in  reference  [21]. 
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1.3  Roadmap  of  the  paper 


For  completeness,  the  perturbation  theory  result  for  the  stochastic  cross-section  per  unit 
area  is  outlined  in  Section  2.  It  is  presented  in  a  way  that  puts  this  approximation  into  the 
context  of  the  formalism  throughout  the  rest  of  the  paper.  Perturbation  theory  is  shown  to 
be  an  approximation  to  the  small  slope  approximation  that  can  easily  be  derived  from  it. 
Perturbation  theory  is  much  easier  to  evaluate  than  the  small  slope  approximation,  and  its 
accuracy  often  suffices  away  from  specular. 

The  small  slope  approximation  is  then  examined  in  the  context  of  several  spectra.  It  is 
assumed  that  a  typical  spectrum  goes  to  zero  for  very  large  scales  (i.e.,  as  the 
wavenumber  approaches  zero),  rises  to  some  peak  and  then  falls  off  as  a  power  law.  For 
wavenumbers  corresponding  to  frequencies  below  below  10  kHz,  the  ocean  surface 
spectrum  is  dominated  by  gravity  waves  (to  about  10  times  the  peak  frequency)  and  then 
enters  a  mixed  gravity-capillary  wave  region  that  extends  down  to  cm  scales.  For  a  part 
of  the  mixed  region  that  extends  up  to  wavenumbers  of  at  least  40  m'1  (corresponding 
roughly  to  the  Bragg  scale  for  10  kHz),  the  power-law  tail  inherited  from  gravity  wave 
still  holds,  and  the  model  above  applies  [22,  23,  24].  (At  higher  wavenumbers,  capillary 
waves  first  cause  the  spectrum  to  drop  off  more  slowly  as  a  function  of  the  wavenumber 
and  then  to  fall  off  sharply-e.g.,  see  references  [22,  24].)  The  ocean  bottom  roughness  at 
scales  relevant  to  acoustic  fields  below  40  kHz  results  from  sedimentation  and  erosion, 
and  this  type  of  rough  surface  again  fits  our  model.  (At  smaller  scales  (i.e.,  larger 
wavenumbers  and  frequencies),  the  ocean  bottom  is  characterized  by  the  presence  of 
small  objects  such  as  sea  shells  and  pebbles,  and  such  objects  generate  power-law  spectra 
that  drop  off  more  rapidly  as  a  function  of  wavenumber.  This  leads  to  a  kink  in  the 
spectrum  at  around  40  kHz  [25].)  Thus,  our  assumption  that  the  spectrum  hits  a  peak  and 
then  falls  off  as  a  single  power  law  is  appropriate  for  the  air-sea  interface  at  acoustic 
frequencies  roughly  below  10  kHz  and  for  the  ocean  bottom  at  acoustic  frequencies 
roughly  below  40  kHz. 

On  the  air-sea  interface,  both  the  tail  and  the  peak  can  generally  be  expressed  in  terms  of 
measurable  physical  quantities  such  as  the  wind  speed.  For  various  types  of  ocean 
bottom,  there  are  criteria  for  detennining  the  power  of  the  drop-off  of  the  tail.  For 
example,  for  exposed  rock  the  wavenumber  typically  drops  off  with  a  power  in  the  range 
2. 6-2. 8  [2,  15],  gravel  with  a  power  3.0  [15],  and  undulating  sand  drops  off  with  a  power 
in  the  range  3.0-3.25  [3,  15].  However,  no  general  criteria  have  found  wide  acceptance 
for  determining  the  peak  of  the  ocean-bottom  spectrum.  It  is  therefore  useful  to  begin 
with  an  examination  of  an  approach  for  calculating  the  stochastic  scattering  cross-section 
when  only  the  behavior  of  the  tail  of  the  spectrum  is  known.  Section  3  outlines  such  an 
approach.  As  outlined  in  references  [2,  20],  these  results  have  been  used  to  deduce 
bottom  parameters  such  as  the  sound  speed,  density  and  roughness  spectrum  by  inverting 
backscattering  results  for  the  ocean  bottom.  These  parameters  can  in  turn  be  substituted 
into  the  theory  to  predict  a  wide  variety  of  monostatic  and  even  bistatic  scattering  results. 

Section  4  presents  the  cross-section  for  a  simple  peaked  spectrum  fonned  by  taking  the 
difference  between  two  power-law  spectra  (denoted  below  as  the  difference  spectrum). 
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This  form  is  useful  in  establishing  general  attributes  of  the  cross-section  related  to  a 
spectral  peak,  and  it  may  accurately  describe  the  spectra  found  along  some  ocean  bottoms, 
but  for  other  applications  the  peak  generated  in  this  way  may  not  be  sharp  enough.  This 
peak  is  not  as  sharp  as  some  found  at  the  air-sea  interface,  and  one  of  the  goals  of  this 
paper  will  be  to  compare  this  spectrum  with  spectra  that  have  sharper  peaks,  and  to 
examine  the  implications  for  the  modeling  of  scattering  from  the  air-sea  interface.  In  fact, 
Section  6  will  demonstrate  that  the  flexible  and  mathematically  tractable  difference 
spectrum  provides  an  adequate  description  of  the  air-sea  interface  despite  its 
comparatively  gentle  peak. 

Section  5  introduces  two  spectra  with  sharper  peaks.  These  two  peaked  spectra  are  more 
realistic,  for  example,  in  considering  the  air-sea  interface.  In  the  first  case,  a  cutoff  is 
introduced,  and  the  spectrum  of  Section  4  is  set  equal  to  zero  below  this  cutoff.  In  the 
second  case,  a  spectrum  very  similar  to  the  Pierson-Moskowitz  spectrum  is  considered. 
The  Pierson-Moskowitz  spectrum  is  modified  so  that  the  integral  giving  the 
autocorrelation  function  has  a  closed  fonn  solution. 

Section  6  considers  the  pure  power  law  spectrum,  the  difference  spectrum  and  the 
modified  Pierson-Moskowitz  spectrum.  After  comparing  the  spectra  directly  and  with 
the  true  Pierson-Moskowitz  spectrum,  this  section  illustrates  their  implications  with  a  few 
air-sea  interface  examples.  The  stochastic  scattering  cross-section  generated  by  the  pure 
power  law  spectrum,  the  difference  spectrum  and  the  modified  Pierson-Moskowitz 
spectrum  are  compared  and  contrasted. 

Section  0  examines  the  stochastic  scattering  cross-section  corresponding  to  an  angle- 
dependent  surface  spectrum.  A  specific  angular  dependence  is  considered,  that  proposed 
by  Elfouhaily  et  al.  [22]  for  the  air-sea  interface.  This  angular  dependence  is  based  on 
the  lowest  term  of  a  Fourier  series  expansion  subject  to  general  constraints  that  must  be 
imposed  on  a  steady-state  rough  surface  spectrum.  It  is  therefore  quite  general,  and 
should  apply  to  a  variety  of  angle-dependent  spectra.  The  angle-dependent  spectrum  is 
the  product  of  an  isotropic  component  and  the  angle-dependent  term.  At  first,  the 
isotropic  component  is  kept  general,  and  the  corresponding  scattering  cross-section  per 
unit  area  corresponding  is  derived.  Then  a  spectrum  that  combines  the  angular 
dependence  of  Elfouhaily  et  al.  with  the  isotropic  difference  spectrum  of  Section  4  is 
considered.  Setting  the  Fourier  coefficient  to  be  a  constant,  and  imposing  another  lowest 
order  expansion,  the  result  once  again  reduces  to  a  sum  of  tractable  one-dimensional 
integrations. 

Section  8  summarizes  the  key  findings  of  this  effort  and  Section  9  provides  some 
concluding  thoughts. 


1.4  The  formalism 

This  section  outlines  the  formalism  used  throughout  this  paper.  Subsection  1.4.1 
discusses  the  meaning  of  a  cross-section  per  unit  area  in  the  context  of  a  quasi-planar 
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rough  interface.  Subsection  1.4.2  derives  the  basic  equations  needed  to  evaluate  the 
stochastic  cross-section  per  unit  area.  Subsection  1 .4.3  discusses  the  parameterization  of 
scattering  geometry-i.e.,  incoming  wave  vector  and  observer  coordinates  (or  equivalently, 
outgoing  wave  vector). 

1.4.1  The  meaning  of  the  scattering  amplitude  and  cross-section  for 
a  quasi-planar  surface 

Consider  a  quasi-planar  2-dimensional  rough  surface  defined  by  z  =  h[x )  ,  where 
x  =  (x,y)  is  a  coordinate  in  a  2-dimensional  plane.  By  convention,  vectors  in  the  de¬ 
plane  will  be  underlined  in  this  paper.  An  incoming  plane  wave  of  the  fonn  elk  r 
(r  =  ( x,  z )  with  z  >  0  is  some  arbitrary  coordinate  in  3-dimensional  half  space)  scatters 

from  this  surface.  The  total  field  in  the  half  space  defined  by  z  >  h(x)  is  given  by  the 
sum  of  the  scattered  wave  ¥-  (r )  and  the  incoming  plane  wave  elkr . 

Following  reference  [6],  the  2-dimensional  Fourier  transform  of  the  scattered  wave  is 
given  by 

®A^z)  =  \XVM’z)e‘qXJ^ 

(2yr) 

and 

X¥Ap)=l°^U’z)e-^- 

According  to  the  Rayleigh  hypothesis,  the  scattered  field  can  be  modeled  by  a 
superposition  of  outgoing  plane  waves,  and  so  :[q,z^  must  be  of  the  basic  form 

(i> z)  =  Vi  (cb  ^  where 

+^k2  -  q-q,  k2>q-q 

+i^jq  ■  q-k2 ,  k2  <  q  ■  q 
and 

'F*  (r  )  =  J Vi  (g> ) eq-e,q‘zd2q  =  J (p-k  (q)eiqrd2q  .  (1.1) 

Formally  taking  the  stationary  phase  as  r  — »oo  (  r  =  |r  | ),  we  find  that  the  scattered  field 
takes  the  familiar  fonn 

^i(k)  =  (27riq7)Vi(q)^r,  (1.2) 

where  now  q  =  k0r  and  k0  is  the  wavenumber  in  the  half  space.  From  Eq.  (1.2),  we  can 
read  off  the  two  forms  of  the  scattering  amplitude  /'■  (f)  and  T^q,k'j , 

T(q,k) 

f-k  (7f)  =  =  (27rick ) v-k  (g) • 


7 


Formally,  the  scattering  cross-section  is  given  by 


a  ■ 


l2 

T(q,k ) 

1 

16;r2 

(1-3) 


As  noted  in  reference  [6],  this  quantity  is  ill-defined  for  an  infinite  deterministic  surface, 
but  if  we  consider  a  translationally-invariant  ensemble  (i.e.,  the  autocorrelation  function 
only  depends  on  the  distance  between  two  points  on  the  surface),  the  stochastic  cross- 
section  divided  by  the  area  (o') /a  is  an  observable  physical  quantity.  This  quantity  will 
be  discussed  next  in  the  subsection  1.4.2. 


1.4.2  The  basic  form  of  the  lowest-order  small  slope  approximation 
for  the  scattering  amplitude  and  cross-section 

The  lowest-order  small  slope  scattering  amplitude  corresponding  to  an  observer  located 
in  the  direction  q/\q\  is  in  general  of  the  form 

T^q,k^j  =  J3  —  ^d2x  e~ *e,Qzh ,  (1.4) 

Qz 

where  Q  =  k  -q . 


The  prefactor  [)  is  an  algebraic  expression,  and  it  is  not  difficult  to  evaluate.  All 
information  about  the  boundary  conditions  resides  in  this  prefactor,  but  it  contains  no 
infonnation  about  the  interface  h(x) .  The  precise  fonn  of  the  prefactor  J3  for  an 

acoustic  field  incident  on  Dirichlet,  Neumann  and  2-fluid  boundary  conditions,  and  for  an 
electromagnetic  field  incident  on  a  perfect  conductor  can  be  read  off  from  the  results 
given  in  Appendix  C.l  of  reference  [6].  The  prefactor  for  an  acoustic  field  incident  on  a 
fluid-elastic  solid  interface  can  be  deduced  from  Eq.  (3.23)  in  reference  [14],  and  the 
details  are  provided  in  Appendix  A  .a 

The  integral  in  Eq.  (1.4), 

—  f  d2x  eiQxei&h 

aJ 

is  the  same  for  all  boundary  conditions.  It  contains  all  the  infonnation  in  T(q,k )  about 

the  surface  h(x) ,  and  it  also  depends  on  the  scattering  geometry,  but  it  is  not  influenced 
by  the  nature  of  the  boundary  conditions  or  the  media  involved. 


a  Using  the  techniques  developed  in  references  [6]  and  [14],  the  prefactor  /?  and  consequently  the 
scattering  amplitude  T^q,k^j  always  manifestly  obey  reciprocity  (i.e.,  T (q,k^  =  T (~k ,-qj  or 

equivalently  T  (@in,@out,<Pout )  =  T  (0ovt,0-m,  (pout  )  ).  Some  other  versions  of  the  small  slope  at  fluid- 

fluid  and  fluid-elastic  interfaces  (e.g.,  that  used  in  reference  [13])  do  not  manifestly  obey  reciprocity  since 
they  are  ultimately  derived  using  an  approach  that  invokes  a  non-reciprocal  version  of  perturbation  theory. 
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The  next  step  is  to  obtain  the  stochastic  scattering  cross-section  per  unit  area  for  a 
translationally  invariant  interface.  The  latter  condition  implies  that  the  autocorrelation 
function  between  two  values  of  the  surface  height  depends  only  on  the  vector  from  one  to 
the  other: 


^/?(x)/*(V))  =  • 

Taking  the  magnitude  squared  and  ensemble  averaging  Eq.  (1.3)  and  using  Eq.  (1.4)  to 
substitute  for  T  [q,  k ) ,  we  have 


cr  = 


\q,k) 


_  W2 


\6n2  16  n2Q] 

Next,  use  Eq.  (36)  of  reference  [6]  (good  for  a  translationally  invariant  Gaussian  surface): 


J  d2x  d2x 


\  JQ(x-x') 


e'aW-H(-)))  =  eXp 


K 

2 


where 

We  now  have 


/j(x-x')  s  2[yl(0)-yl(x-x')] . 


(<7>  = 


14 


16  7t2Q2 


J  d2x  d2x'e 


r-x')  , 

J 

(1.5) 

’*')]• 

(1.6) 

(Note  that  for  the  moment,  no  assumption  of  isotropy  is  made  yet.  This  will  come  below 
just  above  Eq.  (1.10).) 

Perform  the  change  of  variables 

v  =  x-x'  ;  w  =  ^-(x  +  x')  =>  ^ d2x  d2x  =  ^ d2v  d2w , 

and  note  that  the  integral  J  d2w  gives  the  area  a  of  the  ensonified  region,  and  we  have 


cr 


_  14 


-Jrf 


e;  ,  . 

2  -y^v)  iQ'V 

v e  2  e-  . 


„  -  •  0-7) 

a  1  6k  Qz 

Generally  speaking,  the  integrand  fails  to  approach  zero  at  infinity  and  it  is  necessary  to 
extract  a  specular  8  -function  from  the  integral  to  obtain  the  following  alternate  fonn  of 
Eq.  (1-7): 


cr 


14 


a  \6K2Qr 
Even  so,  the  integral 


(In)2  8  (Q)  +  \d2v 


f  q2z  ,  ,  qI  , 
~pH 

z  —e 


iQv 


f  o 2  o2 

e  —e 


iQv 


(1.8) 


(1.9) 
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is  difficult  to  evaluate,  and  this  manuscript  primarily  addresses  this  issue  for  various 
classes  of  y>(y). 


In  all  the  cases  considered  below  except  Section  07,  the  function  /?(y)  is  spherically 
symmetric,  and  Eq.  (1.8)  reduces  to 

(»>_  M2 


a  1  6n2Q2 


{In)1  e  2  P(  ' d[Q)  +  2n^dv 


f  o2  o2  A 

e  2  -e  2 


^0  (0V) 
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.  (1.10) 


This  integral  can  be  evaluated  numerically  provided  that  p  ( y )  and  /?( qo)  are  known. 


1.4.3  Conventions  used  to  characterize  ^  and  q  in  terms  of  angles 

Fig.  2  outlines  the  conventions  used  in  this  paper  to  characterize  k  and  q  . 


Z 


Fig.  2  -  An  incoming  plane  wave  is  characterized  by  a  wave  vector  k  that  lies  in  the  X-Z  plane  at  an  angle 
6m  relative  to  the  X  -axis.  An  observer  is  at  solid  angle  q  j\q  \  =  q /k0  and  labeled  by  the  vertical  and 
azimuthal  angles  6  t  and  <p.tM  respectively.  This  notation  follows  reference  [15].  Note  that 
<pmll  =180  deg  corresponds  to  backscattering. 


These  conventions  lead  to  the  formulas 

k  =  k0  (cos  6mx  -  sin  6mz) 

q=k0  (cos 0out  (cos (poutx  +  sin cpoJ)  +  sin 0oJ ) 


(1.11) 
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Recall  that  k0  is  the  wavenumber  in  the  half  space. 


These  definitions  give  us  the  magnitudes  of  the  horizontal  and  vertical  components  of 
Q=k-q : 


Q:=z\k-q)^Qz=-k0  (sin  0m  +  sin  0out ) 

Q  =  {\-zz)-(k-q^  \q\  =  kj cos2  0m  +  cos2  6>out  - 2cos 0m  cos 9mt  cos 


(1-12) 


<Pa 


1.4.4  Validity  criteria  for  the  small  slope  approximation 


As  its  name  suggests,  the  small  slope  approximation  is  valid  when  the  slope  of  the 
surface  given  by  the  function  h(x)  is  in  some  sense  small.  In  this  section,  we  state  this 
condition  in  more  precise  form. 

The  small  slope  condition  is  that  root  mean  square  of  the  slope  (vh2rms  be  less  than  one. 
Let  us  estimate  this  quantity. 


Often,  the  tail  dominates  the  calculation  of  the  slope.  Therefore,  begin  by  considering  a 
typical  power-law  tail  for  a  spectrum: 


S(/)  = 


wzir  ’ 


(1.13) 


where  w2  is  some  known  scalar  with  units  meters  to  the  4th  power  ( m4 ),  the  scaling 
parameter  h0  is  1  meterb,  and  the  power-law  parameter  is  a  number  above  2  (in 


principle,  there  is  no  upper  limit  on  ).  The  notation  of  Eq.  (1.13)  follows  reference 
[15].  Recall  that  the  integral  I  appears  in  expression  (1.10)  for  the  stochastic  cross- 
section.  For  such  a  tail  (valid  roughly  between  some  peak  wavenumber  lp  and  the 

characteristic  scale  for  Bragg  scattering  |(7|  -higher  wavenumbers  are  transparent  to  the 
field),  we  have 


b  The  parameter  hQ  is  optional.  If  it  is  present,  then  W2  has  units  [length]4,  but  if  it  is  absent  it  has  units  of 
length  to  a  fractional  power:  4  —  .  Since  it  can  lead  to  confusion  in  the  conversion  of  units,  its  use  by 

the  underwater  acoustics  community  has  recently  declined. 


11 


V/r 


2 nw2  r\e  Idl  ■  f 


tail  hr02 


2  7TW2 

~w 


r 


U 


F2 

4-72 


/4-/2 

lP 


4 -y2  4-y2 


In 


0 


\pj 


Yi*  4 


r2  =4 


(1.14) 


The  Bragg  scale  will  in  most  cases  dominate  the  average  slope,  but  for  completeness,  we 
keep  the  contribution  from  the  spectral  peak  and  now  proceed  to  estimate  the  contribution 
from  lower  wavenumbers  between  zero  and  lp.  To  do  so,  approximate  the  bottom  part  of 
the  spectrum  by  the  basic  fonn 

S(l_)  =  AF  0  <l<lp, 

where  A,  /?  >  0  .  Enforcing  continuity  at  /  ,  we  have 


A,t  =  I  A=  3 1  1 


and 


Thus, 


f  ,  V* 


s(L)  =  sm{L) 


\Jp  j 


_  w2  1 

K  ir; 


r  , 


\lpj 


o  <  /  <  / 


V  h 


2 

rms 


2;nv2  r', 


’ow /  hl2F2+p  •'° 


Idl -FF  = 

Jo 


2  nw. 


K%"f 


1 A+4 
_> _ 

P  +  4 


2kw, 


K 

Adding  the  two  contributions  gives  us 

,l4-r2 


/4 -72 

_> _ 

P  +  4 


Vh 


2 

rms 


2;rw, 


h, 


72 


<2 


4-r2 


— / 


4-72 


1 


1 


4-/2  A  +  4, 


v  r2  *  4 . 


Typically,  /?  =  1 .  For  y2  =  4  ,  we  have 


Vh 


2 

rms 


2  7TW, 


In 


\Q 


vA 


+  - 


P  +  4 


(1.15) 


(1.16) 


(1.17) 


(1.18) 


Thus,  the  small  slope  approximation  is  valid  roughly  as  long  as  (vh2^  as  defined  by 

equation  (1.1 7)/(  1.1 8)  is  less  than  one.  This  is  a  crude  criterion,  both  because  the 
calculation  of  the  rms  slope  is  approximate  and  more  importantly  since  the  small  slope 
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actually  depends  on  a  generalized  slope0,  but  it  provides  a  good  first  indication  for  the 
region  of  validity  of  the  small  slope  approximation.  If  this  quantity  is  close  to  but  less 
than  oned,  then  the  higher  order  small-slope  approximation  should  be  considered  [8,  9, 
26],  As  noted  by  Ivakin  [25],  for  typical  ocean  bottoms,  slopes  are  generally  “small”  on 
size  scales  roughly  down  to  the  typical  Bragg  scale  for  acoustic  fields  below  40  kHz. 
Above  this  frequency,  other  methods  must  be  used. 


For  the  pure  power  law,  use  results  (1.17)  with  the  tenn  including  /  dropped  along  with 
(1.18). 


For  completeness,  note  that  the  nns  surface  height  calculated  using  the  same  procedure  is 


i  2 

rms 


2  7TW, 


K 


l2-72 

lp 


1 


-  +  - 


1 


y2-2  P  +  2. 


\ 

Q 

J 

Yi 

2-/2 


(1.19) 


This  result  should  be  used  very  cautiously  because  it  is  much  more  sensitive  to  the 
behavior  of  the  spectrum  at  and  below  the  peak  wavenumber.  The  approximations  used 
in  this  region  are  quite  crude.  However,  other  models  for  the  spectral  peak  give 
qualitatively  similar  results.  For  a  cutoff  version  of  spectrum  (1.13)  used  in  reference  [2] 


S(l)  = 


i+ho  my 


r2/2 


the  rms  surface  height  is  given  by 


h  L  =- 


2/nv, 


For  the  difference  spectrum  given  by  Eq.  (4. 1  )/(4.2),  this  result  becomes 

2/rw2(l-o22) 


/f  = 


blK  (y2  -2)_ 


(1.20) 


(1.21) 


with  a2  a  free  parameter  smaller  than  one  and  h,  given  by  Eq.  (4.1 1)  in  terms  of  the  peak 


wavenumber  /  ,  y2  and  a2.  Note  that  all  these  versions  of  the  result  indicate  that  /znns , 
w2  and  y2  are  not  independent.  Also  note  that  Eqs.  (1.20)  and  (1.21)  are  valid  for  y2  >  2  , 
and  do  not  depend  on  the  Bragg  scale.  In  fact,  Eq.  (1.19)  only  weakly  depends  on  the 
Bragg  scale  fory2  >  2  .  The  dependence  on  the  Bragg  scale  for  the  slope  was  driven  by 


c  The  true  expansion  parameter  for  the  small  slope  approximation  is  in  fact  a  modified  slope  [6,  8,  9].  This 
quantity  roughly  corresponds  to  the  average  surface  height  variation  over  a  correlation  length  divided  by 
the  correlation  length.  This  is  a  difficult  quantity  to  calculate.  Furthermore,  the  expansion  strictly  speaking 
becomes  convergent  rather  than  asymptotic  only  when  the  scattering  theory  is  properly  renormalized  (e.g., 
R.  Dashen  and  G.  J.  Orris,  “Rough  surfaces  and  the  renormalization  group,”  J.  Math.  Phys.,  31  2352-2360 
(1990)).  An  impeccably  accurate  calculation  of  the  effective  slope  would  have  to  take  this  into  account  as 
well. 

d  Subject  to  the  caveats  in  footnote  c,  when  the  expansion  parameter  is  greater  than  one,  then  the  series 
diverges. 
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the  extra  factor  / 2  in  the  integrand  of  Eq.  (1.16).  The  dependence  on  the  Bragg  scale  in 
Eq.  (1 . 17)  is  only  strong  for  y2<  4 . 


2  1st-order  perturbation  theory 


To  obtain  lst-order  perturbation  theory,  take  the  expression  for  (cr)/a  given  by  Eq.  (1.8), 

expand  the  integrand  of  I  in  powers  of  the  autocorrelation  function  A(u)  and  drop  all 

specular  8  -functions  that  appear.  Using  the  convention  in  the  definition  of  the  surface 
roughness  spectrum  given  by  Eq.  (42)  of  reference  [6],  namely, 

<21> 

{In) 

we  are  left  with 

(\If)  =  J  dlY.  A(v)  e'0v  =  (27rf  S(C?)  • 

This  gives  us  the  first-order  perturbation  result  for  the  scattering  cross-section  per  unit 
area 


a  4 


(2.2) 


Note  that  the  prefactor  \(3\  is  similar  to,  but  not  identical  to  that  in  Eq.  (17)  of  reference 
[15].  The  prefactor  |/?|~  manifestly  exhibits  symmetry  between  k  and  -q  . 


The  spectra  considered  in  Sections  3  to  6  are  direction-independent:  )  ■  In 

such  cases,  we  can  use  Eqs.  (1.1 1)  to  evaluate  | Q  : 

\Q\  =  \k-q\  =  kJ cos2  0m  +  cos2  0oM  -  2  cos  0iB  cos 0oul  cos  (poul  . 


3  Cross-section  corresponding  to  a  power-law 
spectrum  with  an  indeterminate  peak 


This  section  evaluates  the  integral  /  defined  in  Eq.  ( 1 .9)  for  a  power-law  spectrum  given 
by  equation  (1.13). 


In  reference  [6],  a  spectrum  of  this  fonn  was  used  to  characterize  the  air-sea  interface.  In 
this  case  (interjecting  a  factor  h^2  to  bring  the  notation  in  line  with  that  used  elsewhere  in 
this  paper) 
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where  «T0BA  is  the  coefficient  of  the  Toba  spectrum,  u  is  the  wind  speed  and  g 
represents  the  constant  of  gravitational  acceleration.  In  reference  [6]  where  =3.5  ,  the 
parameter  aT0BA  was  dimensionless,  but  more  generally  aT0BA  carries  units  of  length  to 
the  power  3.5  -y2  (if  h0  =  1  m  ,  then  it  can  be  absorbed  into  aT0BA  to  give  it  a  power  3.5). 

Since  no  assumptions  are  made  about  the  spectral  peak,  the  results  in  this  section  are 
particularly  useful  in  problems  such  as  scattering  from  the  ocean  bottom,  where  the 
nature  of  the  spectral  peak  is  site  specific  and  rarely  well-known.  References  [2,  17,  19] 
provide  further  information  concerning  the  background  and  application  of  the  expression 
for  the  stochastic  cross-section  derived  here  in  Section  3. 

To  evaluate  expression  (1.9)  for  /  ,  we  will  need  to  evaluate  the  autocorrelation  function 
A (y) ,  which  also  happens  to  be  the  Fourier  transform  of  spectrum  (1.13).  However,  a 
pure  power-law  spectrum  such  as  given  in  (1.13)  has  an  infrared  divergence  (i.e.,  as 
\k\  — >  0  ,  the  integrand  has  a  non-integrable  singularity),  and  so  we  will  have  to 
temporarily  introduce  some  sort  of  spectral  cutoff.  Once  we  have  obtained  an  expression 
for  the  autocorrelation  function  A(v)  =  d  (|y|) ,  we  can  expand  it  in  |y|  and  keep  only  the 

first  tenn.  This  results  in  a  power  law  expression  for  ^4 (| v|) .  It  turns  out  that  this  power- 

law  expression  for  the  autocorrelation  function  corresponds  to  the  only  part  of  the 
spectrum  that  we  knew  in  the  first  place,  the  power  law  tail  given  in  Eq.  (1.13). 
Furthennore,  for  a  power-law  autocorrelation  function,  p(oo)— »oo  (recall  that  Eq.  (1.6) 

gives  us  /?(oo)  =  2[4(oo)-yl(0)]),  and  the  second  tenn  in  the  integrand  of  I  goes  to 
zero  (i.e.,  the  fonn  of  the  integral  given  in  Eq.  (1.7)  is  valid). 

The  integral  that  is  left  is  still  nontrivial.  Its  integrand  is  characterized  by  rapid 
oscillations  bounded  by  a  slowly  decreasing  envelope.  Two  series  expansions  with 
overlapping  domains  of  validity  will  be  used  to  characterize  the  integral. 


3.1  Obtaining  ,i(  y) 


Consider  the  integral  that  results  when  we  take  a  Fourier  transfonnation  to  obtain  the 
autocorrelation  function  A  (y)  from  the  now  cutoff  power  law 


Scutoff  ( /_)  = - 3. - - 

(1+WH)  )■ 


(3.1) 


We  have 


A{v)  =  |  d2l_  S(l_)e,v  =  2n\  dl  lS(l)J0  (, Iv ), 

0 


(3.2) 
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where  v  =  |v|  and  /  =  (/_).  Substituting  Scutott  {]_)  for  S(l_ )  in  (3.2),  we  have 

J0  (Iv)ldl 


A(v)  =  2  7TW2j 


(l 


(3.3) 


Now,  evaluate  (3.3)  using  Eq.  (20)  from  Section  8.5  of  reference  [27]: 


J/M  Jo  (xy)  4*ydx= 


(3.4) 


with 


f{x)  =  xv+A (V  +  a2)  M  . 

In  this  equation,  set  v  =  0  ,  //  =  ;/2/2  -  1  .  and  make  the  substitutions  x  —>  l ,  y —>  v  to 
obtain 

r'A  f  f 


/(')=■ 


r 

r/> 

2\ 

an 

1+1 

V 

\a) 

J 

^2/ 


and 


J0  (Iv)ldl 


v^K^Jav) 


f 

rn 

2\ 

1+ 

V 

J 

22~V  ^2 


fy^\ 

V2J 


(3.5) 


According  to  reference  [27],  Eq.  (3.4)  is  valid  for 

-1  <  v  <  2/j  +  %=>  -1  <0  <  y2  -  y. 

This  condition  is  satisfied  for  physically  relevant  y2>2,  and  the  additional  condition 
a  >  0  will  also  hold  when  we  set  a  =  \/h0  . 

Now,  in  Eq.  (3.5)  pull  out  a  factor  of  a 2  and  use  the  equality  K_v  (z)  =  Kv  (z)  (e.g.,  see 
Eq.  (9.6.6)  of  reference  [28])  and  set  a  =  l/h0  .  Substitute  the  result  into  Eq.  (3.3)  to  get 


a(*)=2«4  J-AF)ML  . 

"(1+(v)2r 


2nw^ 


(  \ 

2 

(  \ 

V 

K 

V 

J 

*-l 

2 

h 

\noJ 

^  (yY 


(3.6) 


A?2 2  r 


V  £  ) 


Next,  use  Eq.  (9.6.9)  from  reference  [28] 


KAz)~\T{v){lzY  asz 


0;  n>0 


to  get 
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and  from  the  definition  of  p(v_ )  =  2[^4(0)-^(v)] ,  we  have 


and  since  Kv  (z)  goes  to  zero  exponentially  as  z  — »  oo  , 

p(°o)  =  2^4(0)  =  —3  — - • 

n0  '2  _ ] 

2 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


3.2  Expanding  to  obtain  a  power  law  autocorrelation  function 

■Mv) 

Following  Section  IIIB  of  reference  [6],  we  next  replace  the  modified  Bessel  function 

K  (v/h0)  with  its  first-order  expansion  to  obtain  a  power  law  for  A (y) .  This  is 

T  1 
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reasonable  because  the  function  Ky  ( v/h0 )  appears  in  the  exponential  in  the  integrand 

T_1 

of  Eq.  (1.9),  and  so  only  small  values  of  v/h0  contribute  to  the  final  result.  Furthermore, 
as  discussed  above  Eq.  (48)  in  reference  [6],  the  first-order  term  in  the  expansion  of 
^4(y)  essentially  corresponds  to  the  power  law  tail  of  the  surface  roughness  spectrum, 

which  is  assumed  here  in  Section  3  to  be  all  we  know  about  the  spectrum.  Note  that  if 
the  power-law  tail  is  extended  arbitrarily  close  to  k  =  0  (i.e.,  the  large-scale  limit),  then 
there  is  no  specular  component  to  the  interface,  and  the  specular  8  -function  vanishes. 


To  expand  Ky  ( vj h0 ) ,  first  use  Eq.  (9.6.2)  of  reference  [28]  to  recast  it  in  terms  of 


Bessel  functions  of  the  second  kind  (now  with  v/  =  iy2-1): 

kAz)=«la±im, 

2  sin  vn 

and  then  expand  for  small  z  using  Eq.  (9.6.7)  of  that  reference: 


(3.11) 


n\T(v  +  n  +  \) 


This  gives  us 


I  (z)  =  f-l  1  +  0(z2)  ;  Zvl(z)  =  - - +  o(z2(v+1)). 

‘vV  ’  {2)  r(i-v)  v  ’  ’  2  r(i+v')  '  > 


2V  r(i  +  v/) 


Note  that  in  our  case,  2[v  + 1)  =  y2  >  2  ,  so  the  next  order  is  C)(z2  j .  Now,  substitute  into 
Eq.  (3.11)  and  use  Eq.  (6.1.17)  of  reference  [28]: 

r(V)r(i-v)  =  — — r(u)sin^v' = 


sin^u 


to  get6 


'■Kv{z)UT{v)T-1- 


nz 


Jlv 


2sin[i/^]2'T(u  + 1) 


r(i-q 

o(r). 


+ 


Substituting  into  Eq.  (3.9) 
p{v)  =  2A(0) 

=  A(0)- 


1- 


2v 


rr(v+i) 

vT(v) 


r(v/)21"1- 


nz 


2sin[v'^]2'T(v  + 1) 


nvz 


2lv  ‘r2  (v  + 1) sin [vtt] 

and  then  recalling  v  =  \y2  -l,  z  =  v/h0  and  the  expression  for  ^4(0)  from  Eq.  (3.10),  we 
get 


cNote  that  Eq.  (9.6.9)  from  reference  [28]  would  only  yield  the  first  term,  and  so  we  needed  to  use  this 
longer  derivation  starting  with  Eq.  (9.6.2)  of  the  reference. 
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p(y)  =  p{y)  =  - 


h22r^T 


3r2  — 

l  2 


-1  \n 


(3.12) 


Note  that  in  this  case  the  integrand  goes  to  zero  at  infinity,  and  so  when  evaluating  the 
stochastic  cross-section  there  is  no  need  to  separate  out  the  specular  8  -function  as  in  Eq. 
(1.8).  We  can  substitute  this  expression  for  y>(y)  directly  into  the  expression  for  the 
cross-section  (1.7).  Performing  the  angular  integration  first,  we  have 


8  nQ]  { 


-fpfi) 

o  z 


■  (|fi|v)  = 


UJ 

\ydy 


And  then  substituting  (3.12),  we  have 

—  =  — ^  .2  f  ydy  exp[-«/2^2]j0(y) 
a  8 nQs  Q\  o 


fpm)J0(y). 


(3.13) 


where 


2  /\2 
n  w2Q: 


qT^1  hr02  2r^2T2  ^ 

~  2 


--1  \n 


(3.14) 


It  is  important  to  keep  in  mind  that  a 
approach  breaks  down. 


■oo  as  y2  — >  4  (or  any  even  number),  and  this 


Also  note  that  the  form  of  Eq.  (3.13)  is  very  similar  to  Eq.  (11)  of  reference  [15]. 
However,  the  coefficient  \/3\  is  not  the  same  as  the  coefficient  i?  A  found  in  that 

reference.  The  differences  between  the  equations  stems  from  the  fact  that  the  equation  in 
the  reference  comes  from  the  Kirchhoff  approximation,  while  (3.13)  comes  from  the 
small  slope  approximation. 


3.3  Evaluating  the  integral  in  Eq.  (3.13) 

The  challenge  remains  to  evaluate  Eq.  (3.13).  In  this  section,  the  integral 

r\  00 

/  =-^-7  =  J ydy  J0(y)exp[-a/2-2]  (3.15) 

Q  o 

is  evaluated  by  replacing  it  with  a  pair  of  power  series  expansions  (recall  that  the  integral 
/  was  originally  defined  in  Eq.  (1.9)).  Subsection  3.3.1  develops  an  expansion  roughly  in 
powers  of  a  that  applies  for  a  <0.15  and  Subsection  3.3.2  develops  an  expansion 
roughly  in  a  x  that  applies  for  a  >0.15  .  The  dividing  line  between  the  ranges  of 
applicability  of  the  two  series  is  somewhat  arbitrary,  but  based  on  experience  obtained  for 
the  special  case  considered  in  reference  [6]. 
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It  should  also  be  noted  that  an  analytic  solution  to  equation  (3.15)  exists  for  y2  =3  (see 

Appendix  B  of  reference  [2]).  This  will  no  longer  be  the  case  when  we  consider 
roughness  spectra  whose  power-law  behavior  is  cutoff  by  a  spectral  peak  (Sections  4  and 
5  below). 


3.3.1  The  expansion  good  for  a  <0.15 


This  section  provides  an  ad  hoc  derivation  for  an  expansion  in  a  that  is  based  on 
expanding  the  exponential  in  Eq.  (3.15).  This  approach  will  necessitate  exchanging  a 
series  and  an  integral  despite  the  fact  that  uniform  convergence  criteria  are  not  met.  The 
resultant  series  is  asymptotic.  A  more  rigorous  derivation  would  follow  reference  [6], 
and  employ  a  Watson-Sommerfeld  transformation.  This  approach  would  also  provide  an 
error  estimate.  However,  the  approach  outlined  here  provides  a  correct  generalization  of 
the  results  in  reference  [6]  with  relatively  little  effort,  and  the  error  estimates  will  be 
obtained  using  insights  from  reference  [6]  and  by  benchmarking  against  a  new  approach 
by  Drumheller  and  Gragg  [29]  (for  a  more  detailed  description,  see  reference  [30]). 


To  allow  an  easier  match  to  fonnulas  in  reference  [27],  make  the  variable  change  y  — »  x 
in  Eq.  (3.15).  Expand  the  exponential 


exp 


-ax 


r2~2 


j_  j(-l )"ccnxnn-2n 


72=0 


n\ 


substitute  into  /  and  exchange  summation  and  integration  (see  previous  paragraphs  for 
caveats)  to  obtain 

*  (-1  fa"" 


i  =i; 


n= 0 


n\ 


|  dx  yfx  . 


,(y2-2)n+} 


J0(x). 


Introduce  an  integrating  factor  e  ax  into  /  and  let  a  — »  0  : 


1  n  =  lim 


a— >0 


j dx  eax \[x  x^1  2^"+^ J0 (x) 


Now,  evaluate  the  integral  using  Eq.  (6)  in  Section  8.6  (p.  29)  of  reference  [27]  with  the 
following  substitutions: 

y  =  1  n  =  n(j2-  2) +  2  v  =  0. 

Also  use  the  Lagrange  polynomial  equality  E°  (x)  =  Pn  (x) .  The  result  is 

f  \ 


1  n  =  lim 


(a2  +1) 


-X[”(^2-2)+2] 


r(«(r2-2)+2  )p,[n_iyi 


(<,=  +  !)' 


(3.16) 


:r(n(r2-2)+2)P<n_2)tl(0)  . 


20 


Next,  evaluate  the  Lagrange  polynomial  Pv  (0)  =  Pn^_2y l  (0) .  Begin  with  Eq.  (8. 1 .4)  of 
reference  [28]  (p.  332) 

K  W  = 


V'tih 


1  V  // 

2~2~T 


as  well  as  Eq.  (15.1.1)  (p.  556) 


F 


f  V  jU  \  V  /LI  1  2^ 

V  2  2  ’  2  +  2  2  ’ 2 ’Z 


F 


l+V—» 
2  2 


■  — 2z- 


1  V  ,  v  F  3_  , 

- - ,  — ,  Z 

v2  2  2  2  2  2 


M 


1  v  // 

—  +  —  — 

2  2  2 


r 


2  2 


to  obtain 


r(c)  ^  r(a  +  n)r(b  +  n)  zn  ,  x 

F (a,b,c;z)  =  .  }  ’  \ f  —  =>  F(a,b,c;0)  =  1 

V  7  r(a)r(Z>)S  r(c  +  n)  n\  y  ’ 

p,P)=K  w= 


fl  vVl 

r 

1  H - 

U  2  J  1 

l  2  J 

Next,  use  Eq.  (6.1.17)  of  reference  [28]:  r(z)r(l-z)  =  ^csc(^z)  with  z  =  . 
Recalling  that  in  Eq.  (3 . 1 6),  v  =  n  ( y 2  -  2 )  + 1 ,  this  leaves  us  with 

1  +  -2{n-2))„ 


r 


H(y2^2)+1 


(»)  = 


r 


V  v 


[(r2-2)j 


and 


r 


/  ,  =  — 


y 

2 

V  L 


l  +  ffo-2) 


r 


jy 


1+>-2>) 

- sin 


3  n  / 

2  +  2<ri“' 

Now,  use  Eq.  (6.1.18)  from  reference  [28]: 

r(2*)=7^22'"KrWr(*+i) 

with  x  =  l  +  n(f2 -2)/2  to  get 


n 


( n  \ 

fL-2) 

2  y 


y 

2 

V  L 


l  +  fL-2) 


_  1  o["(y2-2)+i] 


jy 


r 


A  f 


3  n 


iNL-2)  r  y+y(r2-2) 


y  v 


2  2 


Thus, 


/.~2 

n 


,[«(r2-2)] 


i+yL-2) 


sin 


n 


fn  \ 

>-2) 


and 
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I  -  f.rrfi  ./J.i  ) exp [-«r'  J]-  V 


_v(-iyv 


n= 0 


ft! 


2^(-iy v 

n\ 


(  r 


v- 


“|2 


l  +  H/,-2) 

V  ^ 


sin 


n 


-(r2  -2) 


j[«(?'2-2)] 


and  finally 


-  p2 

^(-l)"«” ' 

r 

(  n  V 

l  +  >2-2) 

-2 

sin 

n\ 

(„  V 

>2-  2) 

1  X 

2[«(y2-2)] 

!  n= 1  n\ 

(  2  J. 

- 

- 

(2  ), 

3  ) 

,(3-17) 


with  «  given  by  Eq.  (3.14).  Note  that  the  n  =  0  term  in  the  summation  would  be  zero, 
and  so  it  has  been  dropped.  Also  note  that  since  the  expansion  is  asymptotic,  we  are  now 
truncating  the  summation  at  some  finite  value.  Based  on  experience  from  reference  [6],  it 
is  generally  safe  to  set  max  «  10 . 


3.3.2  The  expansion  good  for  a  >0.15 


To  obtain  an  expansion  for  /  that  roughly  goes  as  the  reciprocal  of  a  ,  take  Eq.  (3.15) 
and  expand  the  Bessel  function  J0  (y)  in  the  series  in  y  given  by  Eq.  (9.1.12)  of 
reference  [28]: 


*/o(t)  =  1- 


(1!)  (2!)“  tl  («!) 

The  convergence  of  this  integral  is  uniform  (e.g.,  see  proofs  of  Borel’s  theorem). 


We  now  have  the  equality 


i  =  Yj  dy  y2"+1  exP  [~ayr2~2  ]  • 

n= 0  (ft!)  0 


(3.18) 


Now,  make  the  change  of  variable 

z  =  ayYl~2 

to  obtain 


00  1  1  00 
f  /  2n+\  ^  r,-2l  *  f  1 

1  ri-\  i  1 1  nvn  zv  —  1  W-r  -r 

2/1+4-/2 

Yi-2 

vyvF|_  J  - 

0 

a 

~  2n+2~ 
J2~2_ 

<N 

1 

) 

Now,  use  Eq.  (6.1.1)  from  reference  [28] 

oo 

r(x)  =  j dz  zx  le~z  ;  Rex>0 
o 

with  x  =  (in  +  2)/(y,  -  2)  to  evaluate  the  integral  with  respect  to  z  .  Substituting  into  Eq. 
(3.18),  we  get 
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1 
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f  C-T 

f  1  > 

n 

r 

f  2(n  +  \y 

1 

to 

1 

K> 
l _ 

[a2^-2)\ 

6  {n  !)2  1 

U a*‘~2)  J 

{  ri~2  J 

and  substituting  back  into  (3.13)  we  have 

_  w 


2 

1 

r  1 

yH)"| 

f  1  3 

n 

r 

(  2(n  +  iy 

lei 

2 

_Yi -2. 

L  a2/^-2)J 

S(n!)21 

V4a*‘-1') 

l  ^2-2  J 

(3-19) 


with  a  given  by  Eq.  (3.14).  (The  finite  number  of  terms  needed  to  accurately  evaluate 
(3.19)  are  given  in  the  next  section  below.) 


3.4  Summary  for  the  power-law  spectrum  s(l_ )  =  w2/(h0 1 /J) 2 

From  Eqs.  (3.17)  and  (3.19),  we  have 
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where 

Qz=z-(k-q),  Q  =  {\-zz)-(k-q ),  ff(g)  =  .  (3.20) 

[K  |  g ) 

The  terms  of  Q,  and  |(9|  can  be  evaluated  using  Eq.  (1.12).  The  prefactor  /?  is  discussed 
in  Section  1.4.2  and  Appendix  A  .  The  quantities  Nnm  and  «transition  depend  on  y2 ,  and 
are  determined  empirically.  The  results  are  summarized  in  Table  1. 


Table  1  -  Sum-transition  parameter  values  for  Eq.  (3.20) 


Yi 

3.5 

3.6 

3.7 

3.8 

3.9 

^transition 

0.39 

0.25 

0.05 

0.03 

0.025 

N 

max 

6545 

610 

62 

67 

54 

The  branch  good  for  a  <  «transition  is  an  asymptotic  expansion,  and  the  best  choice  for  the 
number  of  terms  is  typically  max  =  10.  The  branch  good  for  a  >  «transition  in  principle 
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always  converges,  but  as  y2  dips  below  3.1,  the  expansion  may  be  so  slowly  converging 
that  it  requires  an  impractical  level  of  precision  and  number  of  terms.  r(x)  is  the 
Gamma  function  defined,  for  example,  in  Section  6  of  reference  [28].  The  scaling 
constant  h0  is  typically  chosen  to  be  1  m.  The  parameter  a  is  given  by  Eq.  (3.14). 

Recall  that  a  — »  oo  when  y2  =  4 ,  and  this  two-series  approach  breaks  down. 

As  discussed  in  references  [29,  30],  the  region  of  validity  of  this  two-series  approach 
turns  out  to  be  limited  to  y-,  >  3 .  The  rational  function  approach  extends  the  region  of 

validity  down  to  y2  =  2.4 .  The  parameter  region  2.4  <  y2  <  3  is  relevant  to  the  study  of 
rocky  ocean  bottoms. 


4  Introducing  a  peaked  spectrum  formed  by  taking  the 
difference  between  two  power  laws 


This  section  evaluates  the  integral  I  defined  in  Eq.  ( 1 .9)  for  a  power-law  spectrum  of  the 
form 
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(4.1) 

(4.2) 


Here,  bx  ,  b2  and  a2  are  chosen  to  fit  some  surface  spectrum.  The  autocorrelation 
function  is  given  by 

A(v)  =  Al(v)-A2(v) 

A{v)  =  \d2lSl{l_)eilv  ;  A,  (v)  =  J  d2l_  S2  {l_)e~iLv  ' 

One  example  is  considered  in  the  discussion  surrounding  Eqs.  (43)  of  reference  [6], 
where  these  parameters  are  adjusted  to  approximate  the  Toba  ocean  surface  spectrum  for 
y2  =  3.5  .  This  section  generates  a  general  procedure  valid  for  arbitrary  y2>  2. 


4.1  Evaluating  the  autocorrelation  function  and  p(v) 

Once  again,  we  need  to  evaluate  Eq.  (1.8),  and  the  first  step  is  to  calculate 
p{y)  =  2(^4(0)-2A(v))  and  p(yo)  =  2^4(0) .  Follow  the  calculation  in  Section  3.1  that 
leads  from  Eq.  (3.2)  to  (3.6)  with  the  substitutions  w2  — >  bx  and  h0  — >  b2h0  to  obtain: 
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Similarly,  the  substitutions  w2  — >  ^  and  /?,,  — >  b2h0/a2  lead  to: 
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Thus, 


^(v)  =  4(v)-4(v) 
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We  also  need  4 ( 0 ) .  To  get  this  quantity,  follow  the  procedure  used  in  Section  3.2  to 
take  the  v  — »  0  limit  of  Eq.  (3.6).  Take  Eq.  (4.5)  and  use  Eq.  (9.6.9)  from  reference  [28] 
to  get  the  first  term  of  the  expansion  for  vj ( b:h0 ) ,  and  then  take  the  limit: 
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(l-a22). 


Finally,  use  the  standard  identity  r(y,/2)  =  ((;/2/2) -  1  j r((y2/2] -  I  )  to  get 

— (1-0?). 
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Substituting  into  the  Eq.  (4.5),  we  get 
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and 
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p(v)  =  2A(0) 


=  2A(0) 


Finally,  we  also  have 
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(4.9) 


4.2  Expressing  the  coefficients  in  terms  of  the  power-law 
notation 

The  spectrum  given  by  Eqs.  (4.1)  and  (4.2)  depends  on  five  variables:  a2,  bt ,  b2,  h0  and 
y2 .  Typically,  the  variables  known  are  the  coefficient  of  the  power  law  tail  w2  (see  Eq. 

(1.13))  as  well  as  h0,  y2 ,  and  the  root  mean-squared  surface  height  £  =  yj(h2}  =  ^4(0)  . 

Additionally,  for  the  sea  surface  in  general  and  sometimes  (but  not  in  general)  for 
specific  locations  on  the  ocean  bottom,  the  wavenumber  location  of  the  spectral  peak  / 
is  known  as  well.  The  spectrum  given  by  Eqs.  (4.1)  and  (4.2)  is  designed  for  this 
contingency.  Then,  only  three  equations  are  needed  in  order  to  express  a2 ,  bx  and  b2  in 
terms  of  w2,  h0,  y2,  l  and  £ .  These  three  equations  are  adapted  from  the  conditions 
listed  in  the  caption  to  Fig.  3  in  reference  [6]. 


4.2.1  Requirement  1 :  Tail  behavior  as  /  -> 


00 


As  /  — »  go  ,  the  spectrum  should  approach  the  asymptotic  form,  Eq.  (1.13): 
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4.2.2  The  spectrum  has  a  maximum  at  l 


The  spectral  peak  occurs  when  the  derivative  is  zero: 
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Or  solving  for  b2 
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(4.11) 


It  should  be  noted  that  for  the  ocean  surface,  a  standard  choice  is  /  =  g/u 2 ,  though  this 

is  a  little  higher  than  the  value  of  the  peak  wavenumber  for  the  Pierson-Moskowitz 
spectrum. 


Note  that  Eq.  (4.11)  can  be  substituted  into  Eq.  (4.10)  so  that  both  bx  and  b2  can  be 
expressed  in  terms  of  a-, .  It  now  remains  to  find  a  criterion  for  determining  a-, . 


4.2.3  Use  the  root  mean-square  surface  height  to  obtain  a2 

A  spectrum  of  the  form  (4.1)/(4.2)  cannot  typically  provide  a  spectrum  that  is  peaked  as 
sharply  as  the  one  implied  by  the  criteria  outlined  following  Eqs.  (43a)-(43c)  in  reference 
[6]  or  as  possessed  by  the  widely-used  Pierson-Moskowitz  spectrum.  To  obtain  the  best 
approximation,  use  a2  close  to  but  less  than  one.  In  reference  [6],  the  choice  a2  =0.9 
was  made  for  the  ocean  surface.  (Section  5  will  introduce  a  method  for  making  the 
spectrum  even  more  sharply  peaked,  and  present  a  technique  for  evaluating  the 
corresponding  cross-section.) 


4.3  Summary  and  numerical  implementation  of  the  difference 
spectrum 

Choose  a2  close  to  but  less  than  one,  then  use  condition  (4.11)  to  obtain  b2 ,  and  finally 
use  condition  (4.10)  to  obtain  bx.  Then,  use  Eq.  (4.6)  to  obtain  ^4(0)  followed  by  Eq. 
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(4.9)  to  get  p{  xj  and  Eq.  (4.8)  to  get  p(y) .  Then,  substitute  the  values  for  /?( oo)  and 
p{v)  in  a  numerical  evaluation  of  (1 . 10)  rewritten  in  the  form 


a 


14 


8  nQ 


j  vdv 


Z  0 


e  1  -e 


-pH 


J, 


(|e|4 


(4.12) 


(Since  the  8  -function  is  zero  for  all  non-specular  geometries  -  i.e.,  \/Q  ^  0  ,  it  is  left  out 
of  Eq.  (4.12).) 


This  integral  may  be  solved  numerically  by  following  the  procedure  developed  in 
reference  [6],  The  integrand  has  a  period  of  n ,  and  so  the  integration  is  performed  for 
the  following  intervals: 

0  ~ >  T0 
To  ->T0+^ 

To  +  ^  — >  y0  +  In 


T0+('Vx-l)^^T0+«max^ 

Experience  gained  in  reference  [6]  indicated  that  nmax  =  40  and  y0  =  were  good 
choices  for  the  air-sea  interface  parameters  used  there.  The  sum  of  the  contributions 
from  these  intervals  can  be  viewed  as  an  alternating  series  of  monotonically  decreasing 
terms.  As  discussed  in  reference  [6]  convergence  of  such  a  series  can  be  accelerated  by 
averaging  adjacent  terms: 

Up  ,  (-l)”"1  un-\  +(-!)”  u„ 

2  2 


More  recent  work  on  the  numerical  evaluation  of  the  integral  in  Eq.  (4.12)  has  modified 
this  picture  [31].  For  many  parameter  regimes  of  interest,  the  behavior  of  this  integral 
differs  profoundly  from  that  encountered  in  reference  [6],  In  particular,  note  that  Eq. 
(4.8)  for  p(v)  suggests  the  change  of  variables  to  x  =  vj{b2h0 )  .  Then,  taken  as  a 


function  of  x ,  the  envelope,  which  can  be  approximated  as1 


2x 

n 

Q 

b2h0  \ 

l  J 

is  primarily  influenced  by  the  parameters  y2  and  a2.  Furthermore,  the  argument  of  the 
Bessel  function  becomes  \Q\b2h0x . 


f  The  square  root  prefactor  comes  from  xJ0  ^  Q  b2h0X^j  ~  ^2 xj 7t  Q  b2h0 


COS  (1^1  ■ b2h0x-a/4)  , 


where  the  cosine  has  been  absorbed  elsewhere. 
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A  typical  curve  comparing  the  envelope  and  the  Bessel  function  is  shown  in  Fig.  3.  It  is 
best  to  perform  a  straight  sum  (i.e.,  without  averaging  adjacent  terms)  out  to  x  =  x  ,  and 
then  sum  pairs  of  alternating  values  for  about  10  additional  terms. 


Fig.  3  -  A  typical  curve  comparing  the  envelope  and  the  Bessel  function  in  the  integrand  of  Eq.  (4.12).  The 
plot  is  of  the  magnitude  of  these  quantities  vs.  the  parameter  X  =  vj ( b2h0  ) .  It  is  best  to  perform  a  straight 

sum  out  to  a  transition  point  X  =  Xp ,  and  then  sum  pairs  of  alternating  values  for  about  10  additional  terms. 
As  y2  varies  from  2  to  4,  the  transition  point  Xp  varies  more  or  less  linearly  from  2.97  to  4.39.  In  this 
curve,  0m  =0oM  =20  deg  ,  (pout  =  180  deg  ,  w=10m/s  ,  /  =  400Hz  ,  £p  =  g/w2  , 

W2  =  \olTOba /(1 00;r ) ]  \ju ~ I g  (with  OCTOBA  =  0.08  )  ,  y2  =  4  and  h0  =  1 .  The  transition  occurs  at 
X p  =  4.39 .  These  parameters  reflect  the  examples  examined  in  Section  6. 


For  a  majority  of  parameter  values,  Q  b2h0  is  not  small,  and  the  number  of  terms  to  x  is 


given  by 


The  total  number  of  terms  is  np  + 10  . 
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In  a  few  cases. 


Q 


b2h0  □ 


1  and  the  envelope  dominates  the  convergence-i.e.,  it  becomes 


very  small  during  the  first  oscillation  of  the  Bessel  function.  In  general,  |()|  is  small  at 
low  frequency,  or  when  0m  and  0out  are  near  90  ,  or  when  0m  and  6out  are  similar  and 
(j)out  is  near  0  .  Likewise,  b2  is  small  for  large  kp ,  which  for  the  sea  surface  occurs  at 
low  wind  speed  U  as,  for  example,  /  □  g/U1  g.  Then,  one  needs  to  integrate  out  to 
where  the  envelope  is  sufficiently  small  (  3tt/4  +  10tt  beyond  xp  is  usually  a  good 
choice).  A  similar  situation  occurs  if  cxp(  -  (J1  p(rsj  )/2)  is  very  small.  Then  the 
envelope  integration  goes  out  to  where  cxp(-02/)(  y  (x)  j/2  )  is  also  very  small. 


For  more  details  on  the  new  numerics,  see  reference  [21]. 


5  Sharpening  the  cutoff 

As  noted  in  Section  4.2.3  above,  a  spectrum  of  the  fonn  (4.1)  (known  as  the  difference 
spectrum)  cannot  be  made  as  peaked  as  some  standard  spectra.  Consider  the  ocean 
surface,  w2  and  y2  are  set  by  the  relatively  easily  observed  tail  of  the  spectrum,  /  is 

widely  acknowledged  to  be  at  least  roughly  in  the  vicinity  of  u2 / g ,  and  the  root  mean 
square  of  the  surface  displacement  is  a  measured  quantity.  These  criteria  are  sufficient  to 
impose  a  peak  on  the  spectrum  that  must  be  sharper  than  is  allowed  by  the  difference 
spectrum.  (As  will  be  seen  in  Section  6,  the  difference  spectrum  in  practice  provides  an 
adequately  accurate  description  of  the  ocean  surface  spectrum.) 

This  section  introduces  two  new  spectra  that  have  sharper  peaks.  Subsection  5.1 
considers  the  spectrum  generated  when  a  lower  cutoff  is  imposed  on  the  spectrum 
considered  above  in  Section  4.  Subsection  5.2  considers  a  spectrum  based  on  the 
Pierson-Moskowitz  spectrum,  but  modified  in  such  a  way  that  the  autocorrelation 
function  can  be  expressed  in  closed  fonn.  The  modified  Pierson-Moskowitz  spectrum 
will  be  explicitly  compared  to  the  true  Pierson-Moskowitz  spectrum  and  to  the  difference 
spectrum  in  Fig.  5  in  Section  6.  The  impact  of  the  modified  Pierson-Moskowitz  spectrum 
and  the  difference  spectrum  on  the  stochastic  cross-section  will  also  be  discussed  in 
Section  6  (along  with  the  impact  of  using  a  pure  power  law  spectrum  and  a  step-function 
cutoff). 


g  In  many  cases,  l  goes  as  U  h,  where  b  is  some  fraction  between  0.5  and  2.  This  result  is  fairly  robust 
as  the  relation  between  /  and  U  is  modified  in  this  way. 
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5.1  A  spectrum  generated  by  sharpening  the  low-frequency 
dropoff  of  the  spectrum  considered  in  Section  4 


5.1.1  A  step-function  low-frequency  cutoff 


A  new  parameter  A  is  introduced  that  sets  a  small  wavenumber  cutoff  for  the  spectrum 
(see  Fig.  1).  The  utility  of  this  spectrum  will  be  addressed  in  Section  6  below.  Here  we 
will  derive  the  corresponding  expression  for  p{y_)  that  must  be  substituted  into  Eq.  (1.8) 

to  evaluate  the  stochastic  cross-section  per  unit  area  ( <j)/a  . 


The  autocorrelation  function  is  given  by 


A{v)  =  \d2l  S(l_)e-iL'- '=  J  d2l  S(l)e-U  v. 

|/J>A 

Since  we  assume  that  the  two-dimensional  spectrum  is  angle-independent,  we  have 
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Note  that  the  cross-hatched  area  in  Fig.  1  is  effectively  removed  from  the  integral. 
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Formally,  the  spectrum  is 
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and  so,  combining  with  Eq.  (5.1),  the  autocorrelation  function  is 

^(v)  =  4(e)  -A2(v),  (5.3) 

with 
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0  0 

The  integrals  from  0  to  infinity  just  reproduce  the  values  for  A1  (y)  and  A2  (y)  from 
Section  4  (from  here  on  renamed  [d,  (y)]  and  [d2  (y)]  ).  Thus,  we  have 
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A(v)=[A(v)] 
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where  [A  ( v  )]  =[A,  ( v  )]  -  [a,  (  v)]  |  is  evaluated  using  Eqs.  (4.3)  and  (4.4).  The 

finite  integrals  will  need  to  be  evaluated  numerically. 

Finally,  we  need  to  evaluate  A  ( 0 ) .  Recall  that  Eq.  (4.6)  provides  a  value  for  the 
autocorrelation  function  in  the  absence  of  a  cutoff:  [.<4(0)]  o .  Introducing  a  non-zero 
cutoff  A ,  we  now  have 
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and  then  substitute  b2  — »  b2  la2  to  get 
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Properly  reduces  to  l—a^  for  A=0 
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Note  that 
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Now,  using  the  definition  of  p{v)  and  Eq.  (5.4): 


(5.6) 


Using  (4.7)  to  evaluate  [^(f)]A_0/[^(C))]A_0  and  Eq.  (5.6)  to  evaluate 

[^(°)]A=o/^(0)’  we  have 


where 
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Properly  reduces  to  X-cij  for  A=0 
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The  finite  integral  from  0  to  A  can  be  evaluated  numerically  using  a  standard  method 
such  as  Gaussian  quadratures. 

Note  that  limfAiv(x)l  approaches  zero  exponentially  and  limJ0(x)— >0,  and  use  Eq. 

x— »oo  L-  v  ' -1  x->oo  v  / 

(5.5)  to  evaluate  A  (0) : 

/?( oo )  =  2a4(0) 


2  7th 


b\K 


Yi 


-1 


(i+(M„a)! 


2  xl1  ^2/2) 


-<z, 


A  /-  x2A<1-^/2)^ 

fb\ 


1+ 


V  a2  J 


(5.9) 


Properly  reduces  to  1-02  f°r  A=0 


In  summary,  for  the  spectrum  with  a  sharp  cutoff,  Eqs.  (5.5),  (5.7)  and  (5.8)  give  p{v_ ) 
and  Eq.  (5.9)  gives  /?(oo) .  This  is  then  substituted  into  Eq.  (1.8)  to  obtain  the  stochastic 
cross-section  per  unit  area. 


5.1.2  An  exponential  low-frequency  cutoff 

Now,  consider  a  spectrum  cutoff  at  low  wavenumbers  by  an  exponential.  In  addition  to 
the  parameter  A  introduced  in  the  last  subsection,  we  need  to  add  two  more  parameters: 

—  c ,  a  measure  of  the  concavity  of  the  relevant  section  of  the  exponential 
curve.  For  0  0 ,  the  spectrum  is  concave  from  the  right.  For  c  »  0  ,  the 
section  of  the  spectrum  between  I  =  b  and  /  =  A  is  very  concave  from  the 
right,  while  for  c  ~  0 ,  the  section  is  nearly  linear.  For  c  <  0  ,  the  spectrum 
is  concave  from  the  left  (see  the  left  bound  of  the  black  area  in  Fig.  1 ). 

-  b  ,  the  x  -intercept  point. 

The  cutoff  parameters  A ,  b  are  both  positive  and  b  <  A  .  The  parameters  A ,  c  and  b 
appear  in  the  spectrum  as  follows: 


34 
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(5.10) 


Si(D 


Recall  that  /  =  /_  .  This  spectrum  is  illustrated  in  Fig.  1 . 

As  before,  the  stochastic  cross-section  per  unit  area  is  calculated  using  the  fonnula  (1.10). 
This  equation  is  integrated  numerically.  To  evaluate  fonnula  (1.10),  we  need  the 
quantities  p[v)  and  p(rfJ)  corresponding  to  spectrum  (5.10).  Since 

p(v)  =  2  ( 4  ( 0 )  -  4  (v)  j ,  it  follows  that  we  must  take  the  Fourier  transform  of  spectrum 

(5.10)  to  find  the  autocorrelation  function  4(v) .  We  need  to  be  able  to  evaluate  this 

quantity  for  arbitrary  v  positive  and  finite,  and  we  also  need  to  determine  its  value  as  we 
take  the  limits  v  — >  0  and  v  — »  go  . 


The  autocorrelation  function  associated  with  this  spectrum  differs  from  the 
autocorrelation  function  generated  in  Subsection  5.1.1  in  that  the  Fourier  transform 
includes  a  contribution  from  the  area  under  the  exponential  (i.e.,  the  area  in  black).  Thus 
Eq.  (5.4)  becomes 
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(5.11) 
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Note  that  from  Eq.  (5.4)  ,  [^4 ( v)]  t  is  the  autocorrelation  function  corresponding  to  the 
step-function  cutoff  discussed  in  Section  5.1.1  above. 
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Next,  we  need  to  evaluate  4(0)  .  Recall  that  Eq.  (5.5)  provides  a  value  for  the 
autocorrelation  function  corresponding  to  the  step-function  cutoff:  [^4 (0)]  (  .  Setting 

v  =  0  in  Eq.  (5.11)  and  noting  that  J0(0)  =  1,  we  obtain  4(0)  corresponding  to  the 
spectrum  given  in  Eq.  (5.10): 
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(5.12) 


Evaluating  the  trivial  integrals  in  Eq.  (5.12)  leads  to 
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Now,  using  the  definition  of  p(v)  as  well  as  Eq.  (5.11)  for  4 (  v)  and  Eq.  (5.13)  for 
4(0),  we  have 
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Equation  (5.7)  gives  us  the  value  of  p  for  the  step  function  p  ( v) ,  which  we  substitute 
into  Eq.  (5.14)  to  obtain: 
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(5.17) 


(The  boxed  quantity  in  Eq.  (5.15)  corresponds  to  Eq.  (5.7),  the  result  /?step(v)  for  the 
difference  spectrum  with  a  step  cutoff.) 

Finally,  recall  that  in  addition  to  p{v),  its  limit  as  v— >oo,  /?(co),  also  appears  in  the 
expression  for  the  scattering  cross-section,  Eq.  (1.10).  To  obtain  this  quantity,  once  again 
note  that  lim  J0  (x)  — »  0,  and  so  A( go)  =  0  .  Use  Eq.  (5.13)  to  evaluate  ^4(0) : 
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We  can  substitute  for  Alcp  (0)  using  Eq.  (5.17): 
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(5.19) 


In  summary,  for  the  spectrum  with  an  exponential  cutoff,  p{\j)  is  given  by  Eq.  (5.15) 

supplemented  with  Eq.  (5.17)  to  give  Alcp  ( 0 )  and  Eq.  (5.16)  to  give  a  and  f)  .  p{ oo)  is 

given  by  Eq.  (5.19).  These  results  are  then  substituted  into  Eq.  (1.10),  which  is  evaluated 
numerically  to  obtain  the  stochastic  cross-section  per  unit  area. 


5.1.3  Discussion:  the  difference  spectra  with  step  and  exponential 
cutoffs 

The  difference  spectrum  with  a  step  function  cutoff  suffers  from  anomalies  and  can  be 
overly  sensitive  to  parameter  values  (e.g.,  see  Fig.  7  and  related  discussion  in  Section  6.2). 
It  should  therefore  be  used  very  carefully  and  sparingly-in  particular  the  cutoff  parameter 
A  should  be  kept  well  below  the  wavenumber  of  the  spectral  peak.  (Recall  that  in  our 
implementation,  A  is  generally  between  the  peak  of  the  difference  spectrum  and  zero.) 
The  difference  spectrum  with  an  exponential  cutoff  is  more  stable,  but  to  avoid  problems 
similar  to  those  encountered  with  the  step  cutoff,  the  exponential  slope  parameter  c 
should  be  kept  small  enough  that  the  exponent  cutoff  to  the  difference  spectrum  does  not 
closely  resemble  the  step-function  cutoff  and  consequently  suffer  from  its  drawbacks.  As 
discussed  in  Section  8,  the  exponential  cutoff  provides  a  peak-sharpening  correction  to 
the  difference  spectrum  in  cases  where  this  effect  needs  to  be  studied  for  non-integer 
values  of  y2 . 

Recall  that  the  exponential  cutoff  is  designed  to  be  inserted  between  the  peak  of  the 
difference  spectrum  and  zero.  Under  these  circumstances,  it  should  be  noted  that  while 
the  exponential  cutoff  provides  a  narrower  peak,  it  does  not  increase  the  height  of  the 
peak  (keeping  the  tail  behavior  fixed).  This  is  the  big  limitation  of  all  cutoff  difference 
spectra.  Some  spectra  for  the  air-sea  interface,  such  as  the  Pierson-Moskowitz  spectrum, 
are  characterized  by  very  large  spectral  values  near  the  peak  relative  to  the  values  along 
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the  power-law  tail.  A  numerically  tractable  spectrum  that  can  in  principle  provide  better 
models  for  this  spectrum  is  examined  next. 


5.2  A  modified  Pierson-Moskowitz  spectrum 


Section  5.2.1  considers  a  modified  Pierson-Moskowitz  spectrum  of  the  form 

(520) 

(V) 

This  differs  from  a  true  Pierson-Moskowitz  spectrum  in  power  of  1//  in  the  exponential. 
This  quantity  would  be  squared  in  a  true  Pierson-Moskowitz  spectrum.  This  spectrum  for 
the  air-sea  interface  is  defined  in  reference  [32]  and  included  here  for  completeness: 


s(L) 


«  J-M/fr4*2)! 

4  xl4 


(5.21) 


with  a  =8.10x10  3  and  /?  =  0.74. 


The  modified  spectrum  is  considered  because  this  spectrum  can  be  Fourier  transformed  to 
generate  a  closed- form  expression  for  the  autocorrelation  function  which  can  in 

turn  be  substituted  into  p(y)  =  2(A(0)- A(y))  in  the  integrand  of  integral  (1.9).  Note 
that  the  spectrum  peak  occurs  at 


a 


For  the  Pierson-Moskowitz  spectrum,  it  turns  out  that  /  » 0.6  g/u1  ,  while  in  the 
discussion  of  the  peaked  spectrum  made  up  of  a  difference  of  power  laws,  we  noted  that 
/  =  g/u2  is  a  standard  choice.  (The  difference  spectrum  is  very  flexible,  and  another 

choice  for  /  could  have  been  imposed  on  this  spectrum  as  well.)  In  the  examples 
considered  in  Section  6,  the  difference  spectra  will  use  the  standard  choice  of  /  =  g/u 2 . 


Section  5.2.2  considers  a  further  modification  to  the  Pierson-Moskowitz  spectrum.  It  is 
characterized  by  a  tail  of  l/ /3  rather  than  the  standard  1  / /4 .  This  spectrum  resembles  a 
spectrum  that  has  been  used  to  model  the  ocean  bottom  [3,  33],  Furthennore,  the  result 
can  in  the  future  be  used  to  benchmark  other  calculations  along  the  lines  of  those 
discussed  in  Sections  5.1  and  6. 
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5.2.1  Standard  l//4  power  law  tail 


The  autocorrelation  function  corresponding  to  spectrum  (5.20)  is  given  by  the  integral 

In  y0(/v) 


V ■/°(/v)= 


2;nw  d 1  r*  „  e 


/?,,  da “ 


J>  -r-'ot'v) . 


The  quantity  in  the  brackets  is  a  Hankel  transform  that  can  once  again  be  evaluated  using 
entry  (15),  Section  8.6  of  reference  [27]: 

Jo  {xy)x/y/dx  =  2y'Ajo  (^y)K  0  (V2qy) » 
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Evaluation  of  the  second  derivative  is  trivial.  Maple8  was  used  to  take  the  derivatives 
and  following  further  simplification,  we  have: 
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(5.22) 


for  the  spectrum  given  in  Eq.  (5.20). 


We  also  have 
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This  gives  us  for  the  “l//4  modified  Pierson-Moskowitz”  spectrum  defined  by  Eq.  (5.20): 
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and 
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The  results  of  Eqs.  (5.25)  and  (5.26)  can  now  be  substituted  into  Eq.  (1.8),  and  the 
integral  I  defined  by  Eq.  ( 1 .9)  can  now  once  again  be  evaluated  numerically. 

Note  that  as  a  — »  0 ,  ^4(v)— >  0  as  \fa 2  (as  does  /?(v)).  Given  that  a  pure  power  law 
like  (1.13)  cut  off  by  a  step  function  has  a  similar  singularity  as  a  — >  0  : 

2;rw,  ,.  r  „  T0  (/v)  _  1 
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this  is  exactly  the  behavior  we  would  expect.  Note  that  Eqs.  (3.7)  and  (3.8)  correspond  to 
a  different  kind  of  low-wavenumber  cutoff  that  does  not  lend  itself  to  a  direct  comparison. 


5.2.2  A  l//3  power-law  tail 

Now,  we  consider  the  spectrum 

m-j. *£■  <527> 

(V) 

Note  that  the  “l//3  modified  Pierson-Moskowitz”  spectrum  defined  by  Eq.  (5.27)  can  be 
used  to  describe  some  rough  ocean  bottoms  or  as  a  one-dimensional  version  of  the  l//4 
case. 
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Spectrum  (5.27)  implies  the  autocorrelation  function 
A (— )  =  “TT^" Jo” dl  ^TJo{,v) 
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It  follows  that 
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This  gives  us  for  the  “  1 / /3  modified  Pierson-Moskowitz”  spectrum  defined  by  Eq.  (5.27): 
p(v)  =  2[zl(0)-zf(y)] 

(5.30) 
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P(  co)  =  2 


(5.31) 
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The  results  of  Eqs.  (5.30)  and  (5.31)  can  now  be  substituted  into  Eq.  (1.8),  and  the 
integral  I  defined  by  Eq.  (1.9)  can  now  once  again  be  evaluated  numerically. 


5.2.3  Issues  related  to  the  numerical  implementation  of  the  modified 
Pierson-Moskowitz  spectrum 

The  numerical  implementation  of  the  modified  Pierson-Moskowitz  spectrum  resembles 
that  for  the  difference  spectrum  described  at  the  end  of  Section  4.3.  The  key  difference 
results  from  the  fact  that  the  parameter  /?(y)  defined  by  equations  (5.25)  and  (5.30) 

involves  Bessel  functions  of  the  first  kind  J0  [y[2av  j  and  J]  [Jlav  'j .  This  modifies  the 

scenario  presented  by  Fig.  3  since  now  the  envelope  itself  oscillates.  Thus,  the  case  of 
|o|/y/?0  needs  to  be  approached  with  increased  care.  One  must  be  sure  to  integrate  out  far 
enough  that  both  oscillations  have  died  down. 


6  A  look  at  the  various  spectra 

This  section  includes  several  figures  illustrating  the  utility  and  acoustic  implications  of 
the  various  spectra. 


6. 1  Comparing  just  the  spectra 

First,  we  include  a  set  of  figures  comparing  the  spectra  themselves.  Examples  of 
parameter  sets  to  be  used  in  the  comparison  follow.  All  expressions  assume  MKS  units 
throughout. 

The  difference  spectrum  (i.e.,  the  spectrum  made  up  of  the  difference  between  two  power 
laws)  is  given  by  Eqs.  (4.1)  and  (4.2).  In  this  section,  the  parameters  in  these  equations 
will  be  assigned  values  appropriate  for  the  air-sea  interface: 
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The  expression  for  bx  and  b2  come  from  Eqs.  (4.10)  and  (4.11)  respectively  with 
/  =  g/u2 ,  while  the  expression  for  w2  comes  from  item  (i)  in  the  caption  of  Fig.  3  in 
reference  [6]. 

The  relevant  numerical  parameters  are: 

«ioba  =  °-08  m3  5-/2 , 

a2  =  0.9, 

y,  is  taken  to  vary  from  2.0  to  4.0, 
g  is  the  gravitational  constant  9.80612  m/s  . 

As  discussed  in  Subsection  4.2.3,  the  value  a2  is  set  close  to  but  less  than  1.  (Note  that 
we  could  also  absorbe  a  factor  of  h^2  into  «Toba  and  have  units  m3'5.) 

Fig.  4  illustrates  difference  spectra  for  the  listed  values  of  y2  (the  power  of  the  tail)  and 
wind  speed.  The  higher  wind  speeds  (corresponding  to  spectral  peak  locations  /  □  1 ) 

considered  in  Fig.  4c  and  d  each  have  a  “cross-over  point”  at  wavenumber  /  =  1 ,  where 
the  order  of  the  largest  to  the  smallest  spectrum  reverses.  This  happens  because  /  =  1  is 

on  the  tail  of  the  spectrum,  so  that  the  spectrum  5(/)|/_i  « w2  j(h0l )/2  ~w2  is 

independent  of  y2  at  this  point  only.  As  the  wind  speed  decreases,  the  spectral  peak 
location  increases  and  the  situation  changes.  For  a  wind  speed  of  5  m/s  (i.e.,  /  ~  0.4), 
then  the  situation  in  Fig.  4b  applies  and  there  is  no  single  cross-over  point.  For  even 
lower  wind  speeds  such  that  /  >  0.4 ,  the  curves  do  not  come  together  at  all  (  Fig.  4a). 
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Fig.  4  -  Difference  spectra  (defined  by  Eqs.  (4.1)  and  (4.2))  with  parameters  appropriate  for  the  air-sea 
interface.  The  tail  of  the  difference  spectrum  goes  as  W2j (/tn/)2  .  y2  steps  from  2.0  to  4.0  in  increments 
of  0.2.  The  wind  speed  varies  from  2  m/s  to  20  m/s. 
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Now  the  difference  spectrum  will  be  compared  to  the  modified  Pierson-Moskowitz 
spectrum  (5.20).  The  new  parameter  in  the  exponent  of  the  modified  Pierson-Moskowitz 
spectrum  is  given  by 

a  =  4lP  =4(&/w2)  (6-2) 

Note  that  g/w2  is  the  standard  peak  wavenumber,  and  this  choice  insures  that  the  peak 
wavenumbers  l  of  the  difference  spectrum  (4.1)/(4.2)  and  the  modified  Pierson- 
Moskowitz  spectrum  (5.20)  coincide.  The  coefficient  w ,  remains  as  above. 

The  difference  spectrum  (4.1)/(4.2),  the  standard  Pierson-Moskowitz  spectrum  (5.21)  and 
the  modified  Pierson-Moskowitz  spectrum  (5.20)  are  illustrated  in  Fig.  5.  Note  that 
compared  to  the  difference  spectrum,  the  modified  Pierson-Moskowitz  spectrum  more 
closely  reproduces  the  sharp  peaks  that  are  accepted  as  characteristic  of  the  air-sea 
interface.  We  have  a  family  of  modified  Pierson-Moskowitz  spectra  that  closely  match 
the  difference  spectra  for  all  wavenumbers  /  above  the  sharp  cutoffs.  In  this  case,  the 
modified  Pierson-Moskowitz  spectra  have  peak  heights  that  are  a  little  lower  than  those 
found  on  the  difference  spectra.  Furthermore,  note  that  the  modified  Pierson-Moskowitz 
spectra  closely  match  the  appropriate  difference  spectra  along  the  tail,  where  the 
dominant  contribution  to  the  scattering  cross-section  occurs. 

It  is  possible  to  use  an  empirical  I"  w,  1  to  modestly  improve  agreement 

r  r  L  z  J Modified  Pierson-Moskowitz  j  r  o 

with  the  difference  spectrum,  but  there  is  no  compelling  physical  reason  to  do  soh.  More 
interestingly,  it  is  possible  to  select  values  for  a  that  produce  a  match  in  peak  height 
between  the  modified  Pierson-Moskowitz  spectrum  on  the  one  hand  and  the  true  Pierson- 
Moskowitz  spectrum  on  the  other,  while  maintaining  the  tail  behavior  of  the  difference 

spectrum.  Specifically,  a  =  0.436(u/15)/' 41p  =  0.436(u/15)/* 4(g/w2)  with  the 
standard  w2  achieves  this  goal,  but  the  modified  Pierson-Moskowitz  spectral  peak  is 
shifted  to  lower  wavenumber  (particularly  at  low  wind  speeds).  It  turns  out  that  such 
modified  Pierson-Moskowitz  spectra  yield  cross-sections  for  one-dimensional  benchmark 
problems  of  the  sort  examined  in  references  [3,  33]  that  are  in  poor  agreement  with  the 
results  generated  by  both  the  difference  spectrum  and  the  standard  Pierson-Moskowitz 
spectrum.  The  latter  two  agree  very  well,  as  does  the  modified  Pierson-Moskowitz 
spectrum  with  the  parameter  a  given  by  (6.2).  It  follows  that  for  scattering  scenarios 
such  as  those  considered  here,  where  the  Bragg  wavenumber  falls  along  the  tail  of  the 
spectrum,  the  peak  location  is  more  important  than  the  peak  magnitude  in  detennining  the 
behavior  of  the  stochastic  cross-section. 


Specifically,  use  equation  (6.2)  combined  with  the  empirical  choice 
tv,  =(l.67-0.0017ir)ri.44-0.29arctan(5fi/-3))lfw,l  to  obtain  closer  overall  agreement 

z  V  /  |_  \\//_|L  z  J  difference  spectrum 

between  the  cross-sections  generated  by  the  difference  and  modified  Pierson-Moskowitz  spectra.  This 
empirical  formula  is  for  /  =  g  /  U2 . 
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WAVENUMBER  I  (rd/m) 


Fig.  5  -  Comparison  of  the  difference  spectra  (4.1)  [solid  curves],  Pierson-Moskowitz  spectra  (5.21) 
[dashed  curves],  and  modified  Pierson-Moskowitz  spectra  (5.20)/(6.2)  [dotted  curves].  The  tails  are 
characterized  by  y1  =  4  power  laws,  and  the  wind  speed  goes  from  2  m/s  to  20  m/s  in  increments  of  3  m/s. 
Note  that  at  low  windspeed,  the  difference  and  modified  Pierson-Moskowitz  spectra  have  tails  that  are 
below  those  for  the  corresponding  Pierson-Moskowitz  spectra.  This  is  an  artifact  of  choice  of  lp  and  the 

definition  of  W2  in  Eq.  (6.1). 
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6.2  Scattering  cross-sections:  Examples 


Fig.  6  provides  examples  comparing  the  cross-sections  associated  with  the  difference 
spectra  and  the  modified  Pierson-Moskowitz  spectra.  The  pure  power-law  spectrum 
given  by  Eq.  (1.13)  is  considered  as  well.  The  stochastic  cross-section  for  the  monostatic 
(i.e.,  backscattered)  case  is  graphed  as  a  function  of  the  grazing  angle,  the  frequency,  and 
the  wind  speed.  We  also  consider  a  bistatic  problem  with  the  incoming  grazing  angle 
fixed  at  20  degrees,  and  vary  the  outgoing  angle.  For  the  pure  power  law,  the  tail  power 
y2  is  chosen  to  be  3.95,  close  to  the  modified  Pierson-Moskowitz  value  of  y2  =4,  but 
not  exactly  equal  to  4.  (Recall  that  the  pure  power  law  is  calculated  using  algorithms 
developed  in  Section  3  that  do  not  work  at  exactly  y2  =  4  (see  the  statement  after  Eq. 
(3.14)),  while  the  modified  Pierson-Moskowitz  spectrum  must  employ  an  integer  power 
law.)  The  difference  spectrum  also  has  a  tail  power  of  y2  =  4 . 

Fig.  6a  fixes  the  frequency  at  400  Hz  and  the  wind  speed  at  10  m/s,  and  graphs  the 
monostatic  cross-section  as  a  function  of  the  grazing  angle.  Fig.  6b  considers  bistatic 
scattering.  The  incoming  grazing  angle  is  fixed  at  20  degrees,  the  frequency  at  400  Hz, 
and  the  wind  speed  at  10  m/s.  The  null  at  160  degrees  is  expected.  This  is  the  angle  for 
specular  reflection,  and  here  coherent  reflected  radiation  dominates.  Fig.  6c  fixes  the 
grazing  angle  at  20  deg  and  the  wind  speed  at  10  m/s.  It  graphs  the  monostatic  cross- 
section  as  a  function  of  the  frequency.  Fig.  6d  fixes  the  frequency  at  400  Hz  and  the 
grazing  angle  at  20  degrees  and  graphs  the  backscattered  cross-section  as  a  function  of 
the  wind  speed. 

The  differences  between  the  results  occur  where  we  expect  them  to  be:  near  specular  for 
Fig.  6a  and  Fig.  6b,  and  at  low  frequencies  for  Fig.  6c.  The  differences  in  Fig.  6d  at  low 
wind  speeds  derive  from  this  particular  choice  of  /  (at  such  wind  speeds,  the  spectral 

peaks  are  at  relatively  high  wavenumbers-e.g.,  for  u  =  1  m/s ,  /  □  10  m'1 ).  Also  note  that 
whereas  a  pure  power  law  becomes  frequency  independent  as  y2  approaches  4,  this  is  not 
true  in  general  for  a  peaked  spectrum. 

Fig.  7  examines  the  stochastic  scattering  cross-section  for  the  difference  spectrum  cut  off 
with  a  step  function  (see  Eq.  (5.2)).  Fig.  7a  graphs  the  bistatic  scattering  cross-section  as 
a  function  of  outgoing  angle  for  a  frequency  of  400  Hz,  a  wind  speed  of  10  m/s  and  an 
incoming  grazing  angle  of  20  deg.  Fig.  7b  magnifies  the  behavior  of  this  graph  in  the 
near-specular  direction  to  highlight  the  unphysical  ringing  associated  with  the  sharp 
cutoff.  This  effect  is  particularly  pronounced  for  a  cutoff  A  wavenumber  near  the  peak 
wavenumber  /  .  Fig.  7c  graphs  backscattering  strength  versus  frequency  for  a  wind 

speed  of  10  m/s  and  an  incoming/outgoing  grazing  angle  of  20  deg.  Note  once  again  the 
the  highest  cutoff  A  produces  pronounced  ringing  at  the  low  frequencies.  Fig.  7d  graphs 
backscattering  strength  versus  wind  speed  for  a  frequency  of  400  Hz  and  an 
incoming/outgoing  grazing  angle  of  20  deg.  In  Fig.  7c  and  Fig.  7d,  note  the  unphysical 
dropoff  in  backscattering  strength  for  high  frequencies  and  wind  speeds,  when  the  cutoff 
is  near  the  spectral  peak. 


48 


Fig.  6  -  Various  curves  comparing  in-plane  scattering  strength  for  three  spectra:  the  pure  power  law  given 
by  Eq.  (1.13)  with  y2  =  3.95  ,  the  difference  spectrum  given  by  Eqs.  (4.1)/(4.2)  with  y2  =  4,  and  the 

modified  Pierson-Moskowitz  (P-M)  spectrum  given  by  Eq.  (5.20)/(6.2)  with  the  tail  fixed  at  y2  =  4 .  Fig. 

6a  illustrates  the  stochastic  backscattered  cross-section  as  a  function  of  the  (incoming)  grazing  angle.  Fig. 
6b  provides  the  stochastic  bistatic  cross-section  as  a  function  of  the  outgoing  angle.  The  incoming  angle  is 
set  to  be  20  degrees,  the  frequency  is  fixed  at  400  Flz,  and  the  wind  speed  at  10  m/s.  Fig.  6c  provides  the 
stochastic  backscattered  cross-section  as  a  function  of  the  frequency.  The  grazing  angle  is  fixed  at  20 
degrees  and  the  wind  speed  at  10  m/s.  Fig.  6d  graphs  the  stochastic  backscattered  cross-section  as  a 
function  of  the  wind  speed.  The  frequency  is  fixed  at  400  Hz  and  the  grazing  angle  at  20  degrees.  (In  all 
plots,  the  outgoing  azimuthal  angle  cpa ,  =  180 deg  •) 


49 


SCATTERED  GRAZING  ANGLE  (deg) 


Fig.  7  -  Scattering  strengths  for  the  difference  spectrum  with  a  step  cutoff  (Eq.  (5.2)).  The  tail  goes  as 
y2  =  4  and  the  step  cutoff  is  at  wavenumber  A  given  by  the  legends.  The  remaining  parameters  are 
defined  by  Eq.  (6.1).  The  legend  in  Fig.  7a  applies  to  Fig.  7b  and  Fig.  7c  as  well.  Fig.  7b  is  a  blowup  of 
Fig.  7a  in  the  vicinity  of  the  specular  direction.  Since  lp=g/u~  and  M  =10  m/s  ,  we  have 

/  =  0.098  m  ,  and  the  parameter  values  for  A  therefore  go  from  the  peak  wavenumber  lp  down  to 
0.00  \I p  .  The  lowest  value  of  A  closely  approximates  the  pure  difference  spectrum.  Note  the  unphysical 
effects  associated  with  a  step  cutoff  A  near  the  peak  wavenumber  lp :  ringing  occurs  for  near-specular 

scattering  and  at  low  frequencies,  and  there  is  an  unphysical  dropoff  in  the  backscatter  at  high  wind  speeds 
and  frequencies. 
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The  anomalous  behavior  exhibited  by  the  cutoff  difference  spectrum  (for  a  cutoff 
wavenumber  A  near  the  peak  wavenumber  /  )  is  associated  with  a  discontinuity.  It  is 

important  to  recall  that  the  difference  spectrum  is  being  cutoff  with  the  step  function  so 
that  the  peak  can  be  narrowed.  For  this  to  be  significant,  it  will  be  necessary  that  A  is 
indeed  near  the  peak.  Thus,  this  implies  a  significant  limitation  of  the  usefulness  of  the 
step-cutoff  difference  spectrum.  The  “softer”  cutoffs  associated  with  the  exponential- 
cutoff  difference  spectrum  (Eq.  (5.10))  and  the  modified  Pierson-Moskowitz  spectrum 
ameliorate  these  problems.  However,  some  care  should  be  used  in  choosing  the 
parameter  giving  the  slope  of  the  exponential  cutoff:  c  (see  Eq.  (5.10)).  An  exponential 
cutoff  with  a  very  large  c  will  approximate  the  step  cutoff  and  consequently  lead  to 
many  of  the  same  problems. 


6.3  Discussion 


The  modified  Pierson-Moskowitz  spectrum  with  the  condition  on  a  (6.2)  leads  to  results 
for  the  cross-section  that  are  similar  to  those  calculated  for  the  corresponding  difference 
spectrum.  The  primary  utility  of  the  modified  Pierson-Moskowitz  is  to  provide  an 
example  of  a  spectrum  that  is  similar  to  the  difference  spectrum,  but  has  a  sharper  peak. 
This  can  be  used  to  verify  that  the  effect  of  peak  sharpness  is  modest  in  such  calculations. 
However,  this  conclusion  will  not  be  true  in  general.  The  reason  that  the  nature  of  the 
peak  had  relatively  little  influence  on  the  stochastic  cross-section  for  the  air-sea  interface 
was  because  these  parameters  generated  a  dependence  mostly  on  the  tail  and  the  region 
immediately  to  the  high-wavenumber  side  of  the  peak.  That  may  no  longer  hold  true  in 
other  examples  such  as  a  rippled  ocean  bottom. 

7  Non-isotropic  spectra 

Detailed  modern  ocean-surface  spectra  such  as  the  Elfouhaily  et  al.  spectrum  [22]  and 
the  Donelan-Pierson  spectrum  [24]  contain  a  dependence  on  the  angle  relative  to  the 
wind  direction1.  Furthermore,  some  swell-related  effects  and  some  ocean  bottoms  can 
generate  spectra  that  depend  on  the  azimuthal  angle.  Therefore,  non-isotropic  spectra 
must  sometimes  be  used.  (For  an  example  of  the  measurement  of  non-isotropic  spectra 


‘As  noted  by  McDaniel  [12],  an  additional  factor  needs  to  be  inserted  on  the  right-hand  side  of  her  Eq. 
(3.9)— i.e.,  Eq.  41  of  [22]-to  obtain  spectra  consistent  with  those  shown  in  the  figures  of  [22].  McDaniel 
[12]  also  provides  an  important  caveat  concerning  air-sea  interface  scattering.  She  notes  that  for  light  wind 
speeds  and  radar  in  the  C-  and  Ku-  bands  (wavelengths  roughly  comparable  to  those  associated  with 
underwater  acoustic  frequencies  above  20  kHz-i.e.,  higher  than  those  being  considered  here),  the  small 
slope  approximation  using  the  directional  spectrum  proposed  by  reference  [22]  gives  poor  agreement  for 
out-of-plane  (i.e.,  azimuthal)  V-V  scattering  from  the  air-sea  interface.  At  these  wavelengths,  capillary 
waves  significantly  affect  the  shape  of  the  air-sea  interface  spectrum,  while  for  acoustic  frequencies  typical 
of  most  long-range  underwater  remote  sensing  applications,  the  surface  spectrum  is  in  a  pure  gravity  region 
or  in  a  mixed  region  that  still  follows  the  power  law  typical  for  gravity  waves  (see  the  second  paragraph  of 
Section  1.3  for  more  on  this).  Finally,  also  note  that  reference  [22]  defines  the  isotropic  spectrum  in  such  a 
way  that  it  must  be  rescaled  by  a  factor  of  1/ 2zr/  to  match  the  conventions  used  here. 
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for  the  air-sea  interface  and  their  use  in  conjuction  with  the  small  slope  approximation, 
see  reference  [10].) 

Obtaining  the  cross-section  per  unit  area  for  a  non-isotropic  spectrum  involves  two  nested 
numerical  integrations.  First,  the  autocorrelation  function  A(y)  is  obtained  using  Eq. 
(3-2)J: 

A(v)  =  \d2ls(l_)eiLv.  (7.1) 


Note  that,  as  indicated  in  Eq.  (3.2),  for  an  omnidirectional  (i.e.,  isotropic)  spectrum,  this 
reduces  to 

oo 

A  (  v )  =  2/rj  dl  lS(l)J0  (/v) .  (7.2) 

o 

Spectra  such  those  developed  by  Donelan-Pierson  and  Elfouhaily  et  al.  have  angular 
dependencies,  and  it  is  desirable  to  find  a  way  to  reduce  the  associated  expression  for  the 
autocorrelation  function  to  a  one-dimensional  form  similar  to  Eq.  (3.2)/(7.2). 

As  will  be  shown  in  Section  7.2  below,  a  spectrum  with  an  angular  dependence  like  that 
of  the  Elfouhaily  et  al.  spectrum  (i.e.,  proportional  to  cos(2^>)  k)  lends  itself  to  a  similar 
closed-fonn  integration.  Subsection  7.2.1  considers  the  general  case,  while  Subsection 
7.2.2  obtains  a  closed-form  solution  for  the  autocorrelation  function  for  a  special  case. 
This  special  case  concerns  a  specific  spectrum  that  combines  the  difference  spectrum 
considered  in  Section  4  with  a  specific  angular  dependence  consistent  with  the  basic  form 
proposed  by  Elfouhaily  et  al.  [22]. 


Once  the  autocorrelation  functions  A  ( 0 )  and  A(y)  are  evaluated,  then  we  have 
/?(y)  =  2(yl(0)-y4(y)).  This  is  once  again  substituted  into  the  fundamental  Eq.  (1.8) 
for  the  cross  section  per  unit  area  (reproduced  here): 
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Recall  that  for  an  omnidirectional  spectrum  this  would  reduce  to  Eq.  (1.10): 


J  It  is  important  to  note  that  Eq.  (3.2)  implicitly  assumes  that  S (  /_)  is  symmetric  with  respect  to  reflection 
across  the  origin  (i.e.,  S  (/J  =  S  /_)  ,  or  in  terms  of  the  azimuthal  angle  <p  ,  S  (/ ,(p  +  ^ )  =  S  (/ ,<p)). 
This  is  inherent  in  the  frozen-surface  approximation,  and  it  is  a  necessary  condition  so  that  and 

sU)  can  both  be  real.  To  put  it  another  way,  a  frozen  (time-independent)  surface  must  be  made  of 
standing  waves,  which  have  equal  backward  and  forward  traveling  components.  Time -dependent  theories 
require  formulations  that  are  sensitive  to  differences  under  the  transformation  I_  — >  —l_ ,  but  this  paper  is 
concerned  exclusively  with  the  time-independent  problem. 
k  (p  is  the  azimuthal  angle-i.e.,  the  angular  coordinate  in  the  X-y  plane. 
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(7.4) 
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Section  7.3  explicitly  considers  a  spectrum  with  the  angular  dependence  proposed  by 
Elfouhaily  et  a/.,  and  outlines  how  in  such  cases,  Eq.  (E8)/(7.3)  can  be  reduced  to  a  more 
tractable  fonn.  Assuming  the  coefficient  of  the  angular  dependence  is  small,  we  obtain  a 
result  like  that  encountered  in  Eq.  (1.1 0)/(7.4)  plus  an  additional  term  involving  a  similar 
one-dimensional  integral. 

This  result  is  quite  general.  The  significance  of  the  specific  spectrum  considered  in 
Subsection  7.2.2  is  that  it  produces  a  closed  fonn  expression  for p(v),  thus  avoiding  the 
need  for  nested  numerical  integrations. 

Section  7.4  will  summarize  the  results  for  non-isotropic  spectra. 


7.1  A  modification  of  the  basic  notation 

Before  we  proceed,  it  is  necessary  to  modify  the  notation  established  in  Subsection  1.4.3. 
Now,  the  positive  x  -axis  points  in  some  preferred  direction  set  by  the  rough  surface  itself 
and  azimuthal  angles  must  be  given  relative  to  this  direction.  The  projection  of  the 

incoming  wave  vector  k  onto  the  x-y  plane  must  now  be  associated  with  an  arbitrary 
direction: 


k=k(  cos6>in  (cos  (pmx  +  sin  (pmy )  —  sin  0mz^ ,  (7.5) 

while  q  is  as  before.  Note  that  in  addition  to  outgoing  azimuthal  angle  cpoul ,  we  now 
need  to  specify  an  incoming  azimuthal  angle  (pm  relative  the  preferred  direction  set,  for 
example,  by  swell  on  the  ocean  surface  or  ripples  on  bottom.  This  is  summarized  in  Fig. 
8. 


7.2  The  autocorrelation  function  for  an  angular  dependence 
proportional  to  cos (2 <p) 


Subsection  7.2.1  derives  the  contribution  to  the  autocorrelation  function  associated  with 
an  angular  dependence  of  the  sort  proposed  by  Elfouhaily  et  al.  [22]  in  that  paper’s  Eqs. 
(49)  and  (67).  In  other  words,  the  spectrum  is  of  the  basic  form 


(7.6) 


with 
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(7.7) 


and 


®{l,<Pi)  =  ^[1  +  A(0COS(2^/)]- 

l_  =  /  (cos  cpjX  +  sin  q>ly'). 


The  function  A(/)  is  a  new  parameter  introduced  in  this  equation. 


z 


Fig.  8  -  The  modified  geometry  for  an  angle-dependent  spectrum.  This  differs  from  Fig.  2  in  that  there  is 
an  intrinsic  preferred  direction  distinct  from  the  direction  defined  by  the  incoming  wave  vector  k  .  Most 
typically  for  the  air-sea  interface,  the  preferred  direction  is  the  direction  toward  which  the  local  wind  is 
blowing.  Flowever,  it  could  in  principle  be  another  direction-for  example,  swell  is  associated  with  far¬ 
away  storms  and  does  not  depend  primarily  on  the  local  wind  direction. 


The  now  angle-dependent  functions  p(v_ )  (and  p(oc) )  needed  to  evaluate  the  stochastic 
cross-section  (Eq.  (7.3))  follows  immediately  from  the  autocorrelation  function. 

Subsection  7.2.2  follows  with  a  derivation  of  the  closed-form  solution  corresponding  to 
the  scenario  where  Aisotropic  (/)  (i.e.,  the  isotropic  component  of  the  spectrum)  is  the 

difference  spectrum  (introduced  in  Section  4)  and  A(/)  is  a  constant. 
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7.2.1  Autocorrelation  function:  The  general  case 


For  the  polar  coordinates  given  by  Fig.  8,  we  have  for  the  spectral  wavenumber  l_  and 
the  vector  y  connecting  a  pair  of  points  on  the  x-y  plane, 

l_  =  /(cos^x  +  sin^v)  ;  y  =  v(cos^vx  +  sin^vv). 


From  Eqs.  (7.1)  and  (7.6),  we  have 

co  n 

A(v)  =  \d2l_  S(L)e'L’-  =2n\ldISlmf,(l)\  dV<5>(l, <p,y- 

0  -71 

To  evaluate  the  angular  integration,  we  need 

1_  =  /(cos^x  +  sin^v)  ;  y  =  v(cos^vx  +  sin^vv), 


(7.8) 


and  so 
So,  we  have 


L'  V  =  Iv[cos<p,  cos^  +  sinp,  sin^v  =  /vcos(^  -<pv). 


ilv  cosip 


A(v)  =  2k j / dl Slsotropic (/)  |  dcpO(l,(pl)e"vm<,pl~^) 

0  —n 

iq)  =  (p,~(pv 

00  K~<Pv 

=  2ft j l dl Sisottopic (/)  J  dy<&{l,y  +  (pv)en 

0  ~K—(pv 

Note  that  the  integrand  of  the  angular  integration  0>(l,(p  + (pv)el,vcos*>  is  2ft  -periodic  and 
drop  the  tilde  to  obtain 

oo 

A  ( v)  =  2n\l  dl  S,„„rk  (/)  S^,  (/, «...  v)  (7.9) 


with 


ilv  cos  (p 


Angular  (^v>V)=  {  ^  0>  (/,  <p  +  )  / 

-7T 

1  n 

=  —  J  d<P  [l  + A(/)cos[2(^  +  ^?v)] 


ilv  cos  (p 


Use  the  sum  formula:  cos[2(^  +  ^v)]  =  cos(2^>)cos(2^v)-sin(2^)sin(2^v)  to  get 
(/.p.,  V)  =  2-  f  d<p  cos  (2p, )  f  d<p  cos  (2«.) 

— ;r  —n 

integrand  is  even 


Next,  use  Eq.  (9.1.21)  of  reference  [28]  to  obtain  for  integer  n 
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inJn  (z)  =  —  f  dcp  cos {n(P)eh™(p . 

2  K  J 

-71 

With  n  =  0 ,  n  =  2  and  z  =  Iv  this  gives  us 

^angular  (/,  <PV,  v)  =  J0  (lv) -  A  (/)  COS  {l<py  )J2{lv),  (7.10) 

and  substituting  into  Eq.  (7.9) 

oo  oo 

A(v)  =  2tt\1  dl Sisotropic  (l)J0(lv)  +  cos (2 <pv )-(-2x)jldl Sisotropic  (/) A (/)  J2  (/v) 

-4 - v - ■  0 - . . . . -.(7.11) 

isotropic  (v)  MV) 

=  4sotropic(V)  +  COS(2  <Pv)A{y) 

Note  that  the  isotropic  part  of  the  autocorrelation  function  ^4is  ic  (v)  matches  the  old 
result  given  by  Eq.  (7.2).  A2  ( v)  is  a  new  quantity  associated  with  angular  dependence. 
Since  J2  (0)  =  0  ,  it  follows  that  A2  (0)  =  0 . 

Now, 

p{y)  =  2(A(Q>)  -A(v))  =  pisotropic  (v)  +  cos(2 <pv)p2  (v)  (7.12) 

with 

oo 

Asotropic  (V)  =  2  (isotropic  (0)  “  4sotropic  (V))  =  J 1  dl  isotropic  i1)  f1  "  J0  (/v)) 

0  .  (7.13) 

00 

Pi  (V)  =  2(>^°T-4  (V))  =  -24  (V)  =  4^J  1  dl  isotropic  (O' A(0J2  (/v) 

0 

The  equalities  on  the  far  right-hand  side  of  Eq.  (7.13)  follow  from  (7.11).  Note  that 
Asotropic  (v)  is  identical  to  the  function  p(v)  corresponding  to  a  pure  isotropic  spectrum 

given  by  Sisottopic  (/) . 

For  any  realistic  spectrum,  two  points  infinitely  far  apart  are  uncorrelated  and 
4so tropic  (°°)  =  A2  (qo)  =  0  .  Thus,  p2  (go)  =  0 ,  and 

P^)  =  Asotropic  l00)  =  24sotrop,c  (O)  •  (7-14) 

Equations  (7.13)  and  (7.14)  constitute  the  core  results  of  this  subsection.  The  next 
subsection  considers  these  results  for  the  difference  spectrum  considered  in  Section  4. 

7.2.2  Autocorrelation  function:  The  difference  spectrum 

This  subsection  evaluates  /?(v)  =  pisot  ic  (v)  +  eos(2^v)yC>2  (v)  corresponding  to  the 
difference  spectrum  (4.1)/(4.2): 
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(i+(w)‘ 


(7.15) 


V 


Furthennore,  the  parameter  A(7)  is  assumed  to  be  a  constant  A,  giving  us 

Q>(l><Pi  )  =  ^[1  +  Acos  (2^/)]' 


(7.16) 


Equation  (4.8)  supplemented  by  Eq.  (4.6)  evaluates  the  first  line  of  Eq.  (7.13)  to  give 

Asotropic  (A?) 


— -1  - 

=  2A  (0)  1-  — _ 

^"isotropic  \A/  1  i  _  2 

A  U  a  *  -  ^ 


«2  2  2  E  At  R  2  1 'Wo 


bnhn  J  1  bji, 


2  V  v2  *0 


(7.17) 


4sotapic(0)  a0- 

2 

Equation  (4.9)  in  Section  4.1  evaluates  Eq.  (7.14)  for  the  difference  spectrum: 

H00)  =  Asotropic  l00)  =  24sotroPic  (0)  • 

The  second  line  of  Eq.  (7.13)  is  all  that  is  left.  It  becomes 

A  0)  =  4 xAb,  (|  (v,  y2 ,  l/[h2/i0  ])  -  £  (v,  /2 ,  a2/[h2/i0  ])) , 

with 

MA 


^(v,y2,a)  =  jldl 


This  is  new.  We  will  need  to  evaluate  an  integral  of  the  fonn 

}  J2{lv)ldl 


(7.18) 


(7.19) 


(7.20) 


(7.21) 


with  a  =\/b2h0  and  a  =  a2/b2h0 .  To  do  so,  use  Eq.  (21)  from  Section  8.5  of  reference 
[27]: 
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\f{x)Jv{xy)y[xydx 

0 

2v+xa2^-v+2T{ju  +  \)T{v  +  \) 

r(y2+y2-/i-i) 


+ 


\  2// -71+3 


r(%-y2+ju+2) 


xFi 


v/ _ A, 

2 


ju,v  +  V,y2a2  /A) 

y2,iu+2-y2-y2-y2a2l4) 


with 


/(x)  =  xA  %  (x2  +  a'>'j  . 

In  this  equation,  set  A  =  2,  v  =  2  ,  iu  =  y2j 2-1,  and  make  the  substitutions  x— »/ 
v  — >  v  to  obtain 


/«=■ 


/ 


/2 


r- 


(p+'T  (/w)' 


r 

rn 

a"2 

i+ 

V 

7 

^2/ 


and 


1  ,^rfe/2-a  (  v/4) 

J  16  T(yj2)  121  2/  7  ' 


2 

rn 

2\ 

1  + 

V 

l a) 

J 

.(7.22) 


vr2“V21 


+  2r^X^2/2  +  l)2)  '  |F7?'i/2;r-/2  +  1’ra/2~1;v2g74) 

Note  that  the  required  condition  -Rev<Re/l<2//  +  %=>-2<2<y2  +  ;)2  holds  for 
physically  relevant  y2  >  2 ,  and  the  additional  condition  a  >  0  will  also  hold  when  we  set 
a  =  l/b2/z0  and  a  =  a2/b2/z0 . 

We  can  obtain  an  alternate  form  for  Eq.  (7.22)  by  making  the  substitution 

rfo/2-2)  l 

r(r2/2)  7^-i¥^.2)' 


Inserting  equality  (7.22)  (with  the  above  substitution)  into  Eq.  (7.21)  gives  us 
%(v,y2,a)  =  77  /  \TA~/  \ ‘  1^2 (2;3-y2/2,3;v  a  /4) 


,  (7.23) 


+ V’ y“r(r!/2 +|/)2) ' 'F> 
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which  is  in  turn  substituted  into  (7.20)  to  obtain  p2  (v) .  This  result  is  combined  with 

(7.17)  for  Asotropic  (v)  to  obtain  p(v)  =  pisotmpic  (v)  +  cos(2<pv)  p2  (v)  for  the  difference 

spectrum  (7.15).  Recall  that  /?( oo)  is  given  by  Eq.  (7.19)  (with  (7.18)).  Thus,  we  have 

all  the  statistical  quantities  that  go  into  the  stochastic  scattering  cross-section 
corresponding  to  the  difference  spectrum  (7.15). 


7.3  The  cross-section  for  an  angular  dependence  proportional 

to  cos(2^>) 


Begin  with  integral  Eq.  (1.7): 


cr 


w 


a  \6n2Ql 
Consideration  of  the  specular  term  proportional  to 


nv  e  2 


p(y)  iQ- , 

e- 


\d2v 


~pH 


iQv 


(7.24) 


(7.25) 


is  deferred  to  the  end  of  this  section  because  it  adds  nothing  new  to  the  problem . 


For  the  angular  dependence  proposed  by  Elfouhaily  et  al.  [22],  Eq.  (7.24)  will  require 
evaluation  of  a  one-dimensional  integral  to  obtain  the  autocorrelation  function111,  and 
another  one-dimensional  integral  to  obtain  the  stochastic  cross-section.  In  principle,  the 
cross-section  also  involves  an  additional  one-dimensional  integral  in  the  angular  direction, 
which  can  be  solved  and/or  tabulated".  However,  here  we  will  choose  to  expand  to  first 
order  in  the  angular  dependent  tenn,  and  this  allows  us  to  evaluate  the  angular  integration 
in  closed  fonn,  winding  up  with  an  expression  very  similar  to  Eq.  (1.10)/(7.4). 


Thus,  consider  the  integral 

l\Q)  =  \d1 IIve 


=  vdv  e 


;  /^isotropic  (^) 


J \d(P 


2  P2{v)cos(2<Pv)  JQv 


I  Note  that  from  Eq.  (7.14),  /?(°o)  is  the  same  as  Asotropic  (°°)>  and  so  adding  the  term  proportional  to 

integral  (7.25)  provides  a  standard  specular  8  -function  that  has  no  angular  dependence,  and  subtracting  it 
guarantees  that  the  isotropic  part  of  the  integrand  in  Eq.  (7.24)  properly  goes  to  zero.  Recall  that  since 

p2  (oo)  — >  0 ,  the  angle-dependent  part  automatically  goes  to  zero. 

m  This  can  be  evaluated  in  closed  form  in  certain  cases.  For  example,  recall  that  Subsection  7.2.2  considers 
the  case  where  the  isotropic  part  of  the  spectrum  *S*isot ropic  (/)  is  the  difference  spectrum,  and  obtains  a 
closed-form  expression  for  the  autocorrelation  function. 

II  The  integral  is  some  function  similar  to  a  Bessel  function  or  a  confluent  hypergeometric  function,  but  it  is 

not  common  enough  to  have  been  officially  defined  and  comprehensively  studied  as  these  have. 
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Now,  with  p(v )  =  -(Q^ 1 2} p2  (v) ,  and 

0  =  |g|(cos^Gx  +  sin^e>')  ;  v  =  v(cos^vx  +  sin^,v) , 

we  have 

I2(p,  \q[<Pq)=]  d  cpv  eHv)cos(2n)^|e|vcos(,v-,e)  (7.26) 

-71 

This  integral  is  a  well-behaved  function,  which  could  be  tabulated  numerically. 
Alternately,  it  could  be  studied  systematically,  and  algorithms  developed  to  solve  it  given 
values  of  p ,  ()|  and  cpQ .  However,  here  we  consider  an  alternative.  We  expand  the 

exponential  to  first  order:  =1  +  /?(v)cos(2^>v) .  Now 

i2{p\^(Pq)  =  J  d(Pv  e'q-  +  p(v)  |d^v  eos(2^)e'|-|vcos(pv'Pc). 

-71  -7T 

The  second  integral  becomes 

^  □□□0DHI  7t-(pQ 

Jd (pv  cos(2 ^ )g'fel'cos(Wi ’’o)  =  j  cos^2^  +  ^e)je'^vcos^ 

-*  -x-Vq 


=  jd^cos(2  (^  +  ^e))e'l?|vC0Sf’ 

integrand  is 
2#  periodic 

^  71 

=  J  dtp  cos^2^  +  ^e^cos(zcos^)  +  /J  dtpv  cos(2(^  +  ^G))sin(zcosi^) 

-Tt  -n 

n  n  _____ — ■ — 

=  cos(2^0)  J d tp  cos(2^)cos(zcos^)-  sin(2^e)  ^dtp^jdaJ^^Jco^z  cos  tp) 


+1  cos 


JL  JL  — - - 

(2 <p0)  j  dtp  cos (icp) sin (z cos  tp)-  zsin(2^0)  fj^_^ia^2^)sm(zcos^) 


=  2cos(2^0) Jd^  cos(2^)cos(zcos^)  +  2zcos(2^e) Jd^  cos(2#>)sin(zcos(5) 

0  0 

Use  Eqs.  (3.715-18)  and  (3.715-13)  of  reference  [28]  both  with  n  =  2  to  get 

I  d^cos(2^)cos(zcos^)  =  ^cos  —  J2  (z)  =  -;rX,  (z) 

0  V  2  J 

r  T2^-3 

I  dtp  cos(2<p)sin(zcos^)  =  ;rsin  —  J2(z)  =  0 

0  V  2  J 

and  recalling  z  =  Q^v 

j  d <pv  cos (2^?v ) e'l?lVC0S(^-?,s)  =  -2^-cos(2^e)  J2  (|g|v) . 

-71 

Thus, 
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and 


1 2  ( P,  |  Q\  ,<Pq  )  =  J  d  <PV  e~ 1  -  In p  ( v)  cos  ( 2(pQ )  J2  ( |  £>|  v) 


h{Q)  =  \°vdv  e  ^/w,(v- 

o2 

—  j*  n  7 Asotropic  (v)  ^Q^- 


2  Asotropic  (^) 


v  e 


J  d (pv  e'-'1 -2;rcos(2 -<Pq)\  vdv  e  £ '  77 

-7T 

Q- 


(v)/2(|0|v) 


e--  + 


old  isotropic  result 

Substituting  into  Eq.  (7.24)  gives 

(cr) 


in'^-zos{2(pQ)\\dv  e-e;»™^p2  (v)  J2  (\Q\v) 


a 


_Al_ 

1 6x2Q 


2 

2^1 


p 


2^,  ^  ^2  Asotropic  ( v)  g/Q  — 


old  isotropic  result 


+;z-£>;  cos(2pG)  J°  vdv  e  -  —  p2 


Asotropic  7) 


q\v) 


(7.27) 


Adding  and  subtracting  the  specular  tenn  (7.25),  and  recalling  that  from  Eq.  (7.14), 
p{ oo)  is  the  same  as  /\otrop]C  (go)  ,  we  have 


a 


\p\ 

2 

(2  x)2eQiP's°mAK)S{Q)  +  \d2v 

(  q2  ^ 

g-^Asotropicfi)  _  e'Tfls°'ropiC(00) 

v  "  y 

iOv 

e-~ 

\67t2Q: 

old  isotropic  result 

+nQ]  cos(2 <pQ)\ydv  e^'^v)  p2  {v)J2  (  q\v) 

.  (7.28) 


Note  that  this  consists  of  the  result  for  an  isotropic  spectrum  given  by  Eq.  (1.1 0)/(7.4)  and 
an  additional  tenn  involving  cp0,  the  azimuthal  angle  between  the  prefened  (e.g.,  wind) 

direction  and  the  projection  of  Q  onto  the  x-y  plane: 

Q  =  £>(cos6q  (cos + sin t»q  v)- sin  6*Qzj 

Q  =  ( 1  -zz)-Q  =  Qcos0q  (cos^Qx  +  sin^Qv)  =  |g|(cos^Qx  +  sin^Qj))  • 

'  S  ' 

The  relationship  between  (pm  ,  (ptml  and  (pQ  is  in  general  complicated,  but  for 
backscattering  it  is  cpm  =  cp0  and  <poLil  =  (pm  + 180  deg . 


Also  note  that  result  (7.28)  does  not  make  any  assumptions  about  SisoU opic(/)  or  A  ( /  ) . 
The  significance  of  the  special  case  discussed  in  Subsection  7.2.2  is  that  this  instance 
provides  closed-form  results  for  y0isotropic  (v)  and  p2  (v) . 
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7.4  Summary  for  non-isotropic  spectra 


Section  0  considers  non-isotropic  spectra-i.e.,  those  that  depend  on  both  the  azimuthal 
angle  and  the  magnitude  of  the  wavenumber.  The  basic  fonn  of  the  spectrum  is  the 
product  of  an  isotropic  spectrum  and  a  term  introducing  the  angular  dependence: 

S{L)  =  2*S, 

(The  notation  is  loosely  based  on  reference  [22].)  In  the  most  general  case,  two  nested 
two-dimensional  integrals  are  required:  one  to  obtain  the  autocorrelation  function  d(y) 

from  the  spectrum  S(l_)  and  a  second  to  obtain  the  stochastic  cross-section.  The 
stochastic  cross-section  involves  an  integral  that  contains  the  function 
p(v)  =  2[4(0)-4(v)J  in  the  integrand. 


However,  if  we  assume  a  specific  azimuthal  dependence  for  the  surface  spectrum-that 
suggested  by  reference  [22]-then  the  problem  simplifies.  Then,  the  first  angular 
integration  can  be  solved  in  closed  fonn  (Eq.  (7.10)),  which  in  turn  leads  to  the  two-part 
expression  (7.13)  for  the  isotropic  component  and  the  coefficient  of  the  angle-dependent 
component  of  the  function  /?(y) .  The  second  angular  integration  leads  to  a  function  that 

can  be  tabulated  (Eq.  (7.26)).  Furthennore,  when  we  expand  the  integrand  in  the  second 
angular  integration  (i.e.,  the  integrand  in  Eq.  (7.26))  to  first  order  in  the  angular 
dependence,  this  integration  can  be  evaluated  to  give  us  the  cross-section  in  terms  of  a 
sum  of  one-dimensional  integrals  (Eq.  (7.28)).  In  fact,  Eq.  (7.28)  involves  isotropic 
terms  that  precisely  mirror  Eq.  (1.1 0)/(7.4)  and  a  new  correction  tenn  associated  with  the 
angular  dependence.  The  new  correction  term  involves  a  one-dimensional  infinite 
integral  that  is  quite  similar  to  the  one-dimensional  integral  appearing  in  Eq.  (1.10)/(7.4). 

In  Sections  3  to  6,  we  examined  various  purely  isotropic  spectra  that  were  chosen  to 
produce  a  closed  form  for  the  quantity  p  (  y )  (or  in  the  tenninology  used  here  in  Section 

7,  /:yolropic  (v) ).  Above  in  Section  7.2.2,  one  of  these,  the  difference  spectrum  (first 
introduced  in  Section  4)  is  used  to  construct  an  angle-dependent  spectrum 
S(l)  =  27rSisotmpic(l)0(l,(Pi)  with  5isotropic (/)  given  by  Eq.  (7.15)  and  0(l,<p,)  by  Eq. 
(7.16): 


S(L) 


bi  \ 


(1 +(b2hjff/1 

f 

1  + 

( h  hl^ 

25 

rVi 

V 

l  U2  J 

J 

[1  + A  cos  (2^)]  . 


(7.29) 


Note  that  angle-averaging  this  spectrum  reproduces  the  old  isotropic  difference  spectrum. 
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The  isotropic  part  of  the  cross-section  given  in  Eq.  (7.28)  can  then  be  solved  (as  before) 
by  numerical  integration  with  Asotropic  (v)  given  by  Eq.  (7. 1 7)/(7. 1 8)  and  Asotropic  (°°) 
given  by  Eq.  (7. 19)/(7. 18).  The  new  angle-dependent  term  can  be  solved  by  numerical 
integration  with  p2  (v)  given  by  Eq.  (7.20)/(7.23). 

This  result  can  be  used  to  benchmark  more  general  numerical  calculations  that  by 
necessity  involve  nested  integrations.  Additionally,  noting  that  most  physically  realistic 
spectra  have  power-law  tails,  the  pure  power-law  result  (summarized  in  Section  3.4)  can 
also  be  used  to  benchmark  the  numerical  evaluation  of  (7.3)  in  the  region  where  Bragg 
scattering  falls  well  within  the  region  of  the  power-law  tail  of  the  spectrum.  Furthermore, 
crude  guidance  about  the  behavior  of  the  cross-section  can  be  obtained  when  isotropic 
results  from  Sections  3  to  6  are  applied  to  one-dimensional  slices  of  non-isotropic  spectra. 


8  Summary  and  discussion 

This  paper  develops  spectra  that  allow  tractable  evaluation  of  the  integral  in  the 
stochastic  cross-section  per  unit  area  (o)/a  given  by  Eq.  (1.8),  the  lowest-order  small 
slope  approximation. 

As  a  first  step,  perturbation  theory  is  generated  by  expanding  the  integral  to  first  order  in 
the  roughness  spectrum  (Section  2).  The  resultant  value  for  the  stochastic  cross-section 
per  unit  area  ( cr) /a  is  given  by  Eq.  (2.2),  where  the  integral  is  simply  proportional  to  the 
surface  roughness  spectrum. 

To  evaluate  the  small  slope  expression  that  involves  the  full  integral,  three  isotropic 
rough  surface  scenarios  are  considered,  and  methods  for  evaluating  the  integral  are 
developed  for  each.  As  discussed  in  Section  1.2,  the  spectra  reflect  one  of  two  tradeoffs: 

—  A  tradeoff  between  the  tractability  of  the  numerical  evaluation  of  the 
cross-section  on  the  one  hand  versus  the  ability  of  the  various  spectra  to 
accurately  reflect  as  many  characteristics  of  observed  physical  roughness 
spectra  as  possible. 

—  A  tradeoff  between  competing  physical  approximations. 

The  three  isotropic  scenarios  considered  are  as  follows. 

-  The  first  scenario  (Section  3)  concerns  an  ensemble  of  rough  surfaces, 
where  only  the  power-law  behavior  of  the  tail  is  known.  In  this  case,  the 
integral  reduces  to  a  pair  of  summations  that  apply  in  different  parameter 
regimes.  This  result  is  expressed  by  Eq.  (3.20)  and  the  associated 
commentary  outlining  the  technique  for  evaluating  (o') /a  . 


63 


-  The  second  scenario  is  discussed  in  Section  4.  It  addresses  peaked  spectra 
formed  by  taking  the  difference  between  two  power  laws.  This  is  known 
as  the  difference  spectrum.  In  this  case,  the  integral  reduces  to  a  one¬ 
dimensional  integral  from  zero  to  infinity  with  a  slowly  decaying 
oscillating  integrand.  This  integral  can  be  evaluated  using  Gaussian 
quadratures  over  finite  intervals  each  covering  one  period  as  the  integrand 
oscillates,  and  the  oscillating  terms  can  be  resummed  to  accelerate 
convergence.  The  summary  section  Section  4.3  outlines  the  numerical 
procedure  for  evaluating  {o)/a  for  such  a  peaked  spectrum. 

-  The  third  scenario  considers  spectra  that  have  sharper  peaks  than  can  be 
constructed  by  simply  taking  the  difference  between  two  power-law 
spectra. 


■  The  first  case  is  discussed  in  Section  5.1.  Here  the  spectrum 
constructed  by  taking  the  difference  between  two  power-law 
spectra  is  artificially  cut  off  at  a  low  wavenumber  in  one  of  two 
ways:  using  a  step  function  cutoff  (Subsection  5.1.1)  or  an 
exponential  cutoff  (Subsection  5.1.2)  grafted  onto  the  spectrum 
cutoff  by  a  step-function.  The  integral  evaluation  described  for  the 
second  scenario  just  above  (the  spectrum  formed  by  taking  the 
difference  between  two  power  laws)  is  now  supplemented  by 
additional  finite  integrations  that  can  be  evaluated  using  standard 
technique  such  as  Gaussian  quadratures.  The  actual  procedures  for 
evaluating  the  (cr) /a  for  these  kinds  of  peaked  spectra  with  an 

additional  hard  low-wavenumber  cutoff  are  summarized  in  the  last 
two  sentences  of  Section  5.1.1  and  the  last  three  sentences  of 
Subsection  5.1.2. 

■  The  second  case  is  discussed  in  Section  5.2,  a  spectrum  is  cut  off 
by  multiplication  to  the  exponential  of  the  reciprocal  of  the 
wavenumber.  The  spectrum  discussed  in  Subsection  5.2.1  closely 
resembles  the  Pierson-Moskowitz  spectrum,  while  the  version 
described  in  Subsection  5.2.2  has  a  different  power  law  tail  (l//3 
rather  than  l//4).  As  with  the  difference  spectrum,  the  associated 
autocorrelation  functions  can  be  evaluated  in  closed  form  and  we 
are  once  again  left  with  a  closed- fonn  oscillating  integrand.  The 
approach  previously  used  to  evaluate  the  cross-section  per  unit 
area  {o') /a  for  the  difference  spectrum  can  be  used  again  in  this 

case.  This  procedure  is  outlined  in  the  last  sentences  of 
Subsections  5.2.1  and  of  5.2.2. 

The  step-function  and  exponentially  cutoff  difference  spectra  provide  the  user  community 
with  extra  tools  for  analyzing  the  roles  of  very  sharp  peaks  (as  perhaps  found  along  a 
rippled  ocean  bottom),  but  their  importance  is  generally  secondary  to  that  of  the  three 
spectra  considered  in  more  detail  in  Section  6:  the  pure  power  law,  the  difference 
spectrum  and  the  modified  Pierson-Moskowitz  spectrum. 
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The  power  law  spectrum  without  any  peak  is  primarily  useful  in  cases  where 


—  little  is  known  about  the  nature  of  the  spectral  peak,  while  the  nature  of  the 
tail  is  known, 

-  and  the  Bragg  scale  is  well  into  the  power-law  region  of  the  spectrum,  and 
the  influence  of  the  peak  on  the  scattering  cross-section  is  relatively 
modest. 

This  scenario  applies,  for  example,  to  the  ocean  bottom  in  many  parts  of  the  world.  It 
should  be  noted  that  the  peak  can  be  measured  and  recorded  for  specific  locations.  In  this 
case  the  difference  spectrum  and  possibly  its  “peakier”  variations  will  prove  particularly 
useful,  since  there  are  a  wide  variety  of  y2  ’s  involved.  Application  of  even  modest 
cutoffs  can  sometimes  significantly  alter  the  scattering  (see,  for  example,  Fig.  6). 

Section  6  examined  the  various  spectra  in  the  context  of  acoustic  scattering  (up  to  about 
10  kHz)  from  the  air-sea  interface.  These  examples  relied  on  numerical  implementation 
of  the  theory  by  Joe  Fialkowski  (NRL).  They  indicate  that  the  extra-sharp  cutoff 
introduced  by  the  modified  Pierson-Moskowitz  spectrum  has  a  relatively  modest  effect 
on  the  cross-section  generated  by  the  difference  spectrum.  They  also  indicate  that  the 
pure  power-law  spectrum  can  introduce  some  error,  particularly  at  the  high  grazing 
angles,  near  the  specular  reflection,  and  at  low  frequencies  and  wind  speeds.  The 
difference  spectrum  provides  the  best  balance  between  the  tradeoffs.  The  hard  cutoff  at 
A  <  /  on  the  other  hand  stands  out  as  a  poor  choice. 

Section  0  considers  the  issues  associated  with  evaluation  of  the  stochastic  cross-section 
corresponding  to  angle-dependent  two-dimensional  spectra.  The  angular  dependence 
proposed  by  Elfouhaily  et  al.  [22]  is  examined  in  detail.  The  resultant  cross-section  Eq. 
(7.28)  takes  the  fonn  of  the  isotropic  result  plus  a  correction  associated  with  azimuthal 
dependence.  The  basic  approach  used  to  calculate  the  cross-section  associated  with  the 
isotropic  spectrum  can  be  applied  to  evaluate  the  new  term  associated  with  the  azimuthal 
dependence.  This  is  done  explicity  for  the  spectrum  consisting  of  the  product  of  the 
difference  spectrum  and  the  angular  dependence  proposed  by  Elfouhaily  et  al.  The  result 
is  Eq.  (7.20)/(7.23),  which  provides  the  quantity  needed  so  that  the  new  integral  that 
appears  in  Eq.  (7.28)  reduces  to  the  now  standard  numerical  calculation  of  a  one¬ 
dimensional  infinite  integral  with  an  oscillating  integrand. 


9  Conclusions 

This  work  presents  the  formalism  associated  with  numerical  solutions  of  the  incoherent 
cross-section  per  unit  area  for  quasi-planar  stochastically  rough  surfaces.  The  lowest- 
order  small  slope  approximation  is  widely  used  to  model  the  scattering  cross-section  per 
unit  area  from  such  interfaces.  This  approximation  consists  of  a  prefactor,  which 
contains  information  about  the  scattering  geometry,  the  media  involved  and  the  boundary 
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conditions,  and  an  integral  that  contains  information  about  the  scattering  geometry  and 
the  surface  roughness  spectrum.  The  former  involves  a  closed-form  algebraic  expression, 
but  evaluation  of  the  integral  is  nontrivial.  Some  simplification  of  this  integral  occurs 
when  the  surface  is  assumed  to  be  isotropic,  but  for  this  integral  to  be  truly  tractable,  the 
isotropic  surface  spectrum  must  be  chosen  with  care.  The  choice  primarily  involves  a 
tradeoff  between  tractability  and  the  completeness  of  the  physical  description. 

Generally,  the  most  useful  expressions  for  the  stochastic  scattering  cross-section  come 
from  perturbation  theory  (where  the  integral  disappears  altogether)  and  the  small  slope 
result  associated  with  the  difference  spectrum  (where  the  integral  is  numerically 
tractable).  The  former  is  extremely  tractable:  it  is  trivial  to  calculate  for  any  spectrum. 
However,  it  must  be  used  with  care.  As  noted  in  Section  6  above  and  in  references  [6] 
and  [2],  the  associated  cross-section  for  scattering  from  the  air-sea  interface  and  the 
ocean  bottom  may  in  some  cases  differ  significantly  from  the  correct  value. 

The  small  slope  result  associated  with  the  (isotropic)  difference  spectrum  is  also  very 
flexible:  it  works  for  any  power  y2  >  2 .  The  peak  location  can  be  selected  at  will,  and 
the  spectrum  shape  can  be  adjusted  from  narrow  to  wide.  Physical  spectra  found  in 
underwater  acoustics  applications  typically  obey  2  <  y2  <  4  °.  Even  though  the  modified 

Pierson-Moskowitz  spectrum  (characterized  by  a  very  sharp,  tall  peak)  provided  the  most 
faithful  rendering  of  the  known  characteristics  of  the  air-sea  interface  spectrum,  it  was 
found  that  the  difference  spectrum  provides  perfectly  adequate  values  for  the  scattering 
cross-section.  It  can  be  concluded  that  for  roughly  similar  frequencies  and  spectra  (as 
may  occur  in  acoustic  scattering  from  the  ocean  bottom),  the  difference  spectrum 
provides  the  optimal  tradeoff  between  tractability  and  broad  applicability.  Hence,  it  is  a 
good  default  spectrum  to  use  for  underwater  acoustics  applications  in  the  low-  to  mid- 
frequency  range. 

A  step  function  or  exponential  cutoff  can  be  grafted  onto  the  basic  difference  spectra.  In 
a  typical  case,  this  cutoff  will  be  placed  between  the  peak  of  the  difference  spectrum  and 
zero  wavenumber.  Such  spectra  can  provide  a  platfonn  for  examining  the  effects  of  peak 
narrowness,  but  the  cutoff  procedure  does  not  by  itself  affect  the  height  of  the  peak.  This 
limits  the  utility  of  these  cutoff  difference  spectra,  since  it  is  generally  acknowledged  that 
some  physical  spectra  (e.g.,  the  Pierson-Moskowitz  spectrum  for  the  air-sea  interface)  are 
characterized  by  spectral  values  near  the  peak  that  are  very  large  relative  to  the  values 
along  the  power-law  tail.  It  should  be  further  noted  that  the  step  function  cutoff  generates 
a  very  strong  dependence  on  the  choice  of  parameters,  and  that  therefore  this  spectrum 
should  be  used  very  carefully  if  at  all. 


0  Spectra  with  power-laws  y2  <2  yield  autocorrelation  functions  that  are  extremely  sensitive  to  some 

high-wavenumber  cutoff  imposed  either  by  the  field  (the  renormalization  problem)  or  by  the  characteristics 
of  the  physical  interface  itself.  If  they  are  employed  in  perturbative  scattering  theories  (including  the  small 
slope),  extreme  care  must  be  exercised  to  ensure  that  key  expansion  parameters  are  small.  Spectra  with 
power-laws  y2  >  4  are  perfectly  well-behaved  at  high  wavenumbers  and  should  produce  decent  expansion 

parameters,  but  they  are  quite  sensitive  to  low  wavenumbers.  There  are  no  examples  of  such  spectra  in  this 
paper. 
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The  modified  Pierson-Moskowitz  is  a  very  useful  tool  for  examining  the  effects  of  peaks 
that  are  both  narrower  and  taller  than  is  allowed  by  the  difference  spectrum  (such  as  on 
the  air-sea  interface  and  possibly  a  rippled  ocean  bottom).  While  the  peak  location  can 
be  chosen  at  will,  the  modified  Pierson-Moskowitz  spectrum  only  provides  a  closed  form 
auto-correlation  function  (and  consequently  a  tractable  cross-section)  when  the  tail  goes 
as  an  integral  power  {y2  =3  or  4).  The  width  of  the  peak  follows  once  the  location  of 
the  peak  is  chosen.  For  the  difference  spectrum  the  height  of  the  spectrum  is  fixed  by  the 
behavior  of  the  tail,  the  choice  of  the  peak  wavenumber,  and  the  shape  of  the  spectrum. 
For  the  modified  Pierson-Moskowitz,  the  shape  is  fixed  so  that  the  first  two 
considerations  alone  fix  the  height  of  the  peak.  For  the  modified  Pierson-Moskowitz 
spectrum,  the  height  of  the  peak  can  vary  widely  depending  on  the  value  of  the  peak 
wavenumber.  This  rigid  adherence  to  a  limited  set  of  peak  shapes  and  to  integer  power 
tails  limits  the  utility  of  the  modified  Pierson-Moskowitz  spectrum.  It  can  certainly  be 
productively  used  as  a  benchmark  (here  used  to  verify  the  adequacy  of  the  difference 
spectrum),  and  and  perhaps  for  modeling  the  air-sea  interface  with  a  tail  that  goes  as  the 
integer  4P.  Some  rippled  ocean  bottoms  may  go  as  the  third  power,  and  the  modified 
Pierson-Moskowitz  spectrum  may  also  be  used  to  model  the  scattering  from  such  bottoms. 

The  angular  dependence  proposed  by  Elfouhaily  et  al.  [22]  generates  a  new  tenn  that 
must  be  added  to  the  isotropic  cross-section.  This  new  term  is  very  similar  to  the 
isotropic  result.  The  methods  were  used  to  calculate  the  cross-section  associated  with  an 
isotropic  spectrum  can  be  therefore  be  adapted  to  evaluate  the  new  term  associated  with 
the  azimuthal  dependence. 

The  basic  forms  of  all  these  results  apply  to  both  acoustic  and  electromagnetic  fields  (the 
latter  subject  to  the  caveat  noted  in  footnote  i  on  page  30). 
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p  It  has  long  been  an  unsettled  question  whether  the  tail  characterizing  the  air-sea  interface  decays  as  l  4  or 
whether  it  decreases  as  some  fractional  power. 
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Appendix  A  The  prefactors  corresponding  to  Dirichlet 
and  fluid-elastic  interfaces 

This  appendix  discusses  the  evaluation  of  the  prefactor  (i  found  in  the  small  slope  and 
perturbation  theory  results  for  scattering  from  a  rough  Dirichlet  or  fluid-elastic  solid 
interface.  Although  it  is  associated  with  the  rough  surface  small  slope  scattering 
approximation,  ft  in  general  depends  on  the  quantities  associated  with  a  plane  wave 
scattering  from  a  planar  fluid-solid  interface. 

A.  1  Fluid-solid  boundary  conditions 

The  rough  fluid-elastic  interface  will  be  considered  first.  The  label  I  will  be  connected 
with  the  fluid  side  and  II  will  be  associated  with  the  elastic  solid.  The  incident  and 
scattered  fields  are  evaluated  in  the  fluid.  The  physical  quantities  input  into  ft  are  the 
angles  6m ,  6oat  and  (poM  defined  in  Fig.  2  in  Section  1.4.3,  the  densities  of  the  media  pt 
and  pu ,  the  sound  speed  in  the  fluid  c0 ,  the  Lame  parameters  for  the  solid  Xu  and  pn , 
and  the  frequency  of  the  field  co  .  From  these  parameters,  we  get  the  primary(  p  )-  and 
secondary(  5  )-wave  sound  speed  and  wavenumbers  in  the  solid: 


as  well  as  the  incident  wavenumber  k0  =co/c0. 

The  incident  wavevector  can  be  written  in  terms  of  horizontal  and  vertical  components 
ki  =k  +  kizz  (recall  that  vectors  in  the  horizontal  plane  ( x-y  )  are  underlined),  where 
k  =  k0  costfinx  and  k0z  =  -k0  sin  6m  .  The  sign  in  front  of  sin  6in  reflects  the  conventions 
in  Fig.  2.  The  transmitted  wave  for  the  flat- interface  problem  is  characterized  by  the  p  - 
and  5  -  wavevectors  given  by: 

K=k  +  kJ  kp:=-Jk;-\k\2 

K=k  +  kj  ksz  =  ~^k;-\k\2  •  (A 1 ) 

\kf  =  k20  cos2  em 
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The  prefactor  /?  also  depends  on  the  conjugate  problem.  This  problem  involves  the 
incoming  wavevector  k'  =  -q  =  -k0r  (where  r  is  the  direction  of  the  observer): 

q  =  -k)  +  k{]  sin  8oiaz 


K  =  k'~k0  sin  8outz 


Z=k'+k'pj 


K>=-K 


M 

K  =  k'  +  k'J  k'sz  =  -yjk'2-\k'\2 

K  =  ~K  cos ^out  (cos  ^out*  +  sin  VoJ)  \kf  =  K  cos2  6>out 

The  prefactor  /?  is  deduced  from  Eq.  (3.23)  of  reference  [13]  with  h  — >  z  : 

out  ’  Pi  ’  PlI  ’  K ,  kp  ,  ks  j  — 


(A2) 


4  1-^ 

l  Pi) 


Kqi2 


(ai  +l)(a2  +1) 


+  4 


+8 


Pii 

piK 


( k-q )2 


(ai  +  l)(a2  +  1) 


|  _  Pii 

V  Pi) 

-4k -q 


k-q~,  w  , 

-(ai  +l)(a2  +  1) 

(aj  +l)(a2  +  1) 


(A3) 


+4  k: 


-4k: 


-  +  2£ 


2  Pi  ( al-bl)(a2-b2 ) 
(a1+l)(a2+l)  pu  (a1  +  l)(a2  +  l) 

bpi 


(aj  +l)(a2  +1) 


The  label  “1”  refers  to  quantities  derived  from  the  flat-interface  solution  generated  by  the 
incoming  plane  wave  k  ,  and  the  label  “2”  refers  to  those  same  quantities  evaluated  for 
an  incoming  wavevector  for  the  conjugate  problem  k'  =  -q . 

We  are  now  ready  to  read  off  the  quantities  ax ,  a2,  bx,  b2,  cj  and  g2 .  We  have 


a(d;pnpn;k0,kp,ks) 


pn  k0  sin  0 


Pu_K 
Pi  K 


4|  if 

k4 


4k 2  cos2  8^jk2  -  kl  cos2  6  ^k2  -  k 2  cos2  6  -  ^ jk 2  -  k( 


-  k2  cos2  6 


Pi  yjk2  -  kl  cos2  6 

a\  ~  aiPm  •>  Pn  Pu  ’  ’  kp  ■>  ks )  a2  —  ci  ($out ,  p: ,  pn ,  k0 ,  kp ,  ks  ^ . 


and 

Similarly 

and 


(A4) 


Pu 

K 

r 

i  _ 

2  k2 

_  Pu_ 

k0  sin  6 

Pi 
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^  KJ 
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s Jk 2  -  k2  cos2  8  _ 
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P  k ( ^in  ’ Pi  ’ P u ’ K ’ kp , ks  j  b2  b[6out, pn pn,ka,kp,ks^ .  (A5) 

We  also  have 


and 

Zl=4(0n>K  ’kp’ks)  %2=Z(0Out’ko’kP’ks)-  (A6) 

We  can  also  read  off 

K  =  ~K sin  °m  ;  q0z  =  -^0' =^osin^ut>  (A7) 

and  finally  use 

k=kQcosOmx  ;  -q_  =  k'  =  -k()  cos0oul  (cos^oulT  +  sin ^outy) 

to  obtain 

k-q  =  k-(\-zz)-q  =  k1-cosdmcos 0out  cos ^out .  (A8) 

Note  that  for  pure  (i.e.,  monostatic)  backscatter,  we  have  6m  =  Oom,cpom  =  180° . 

Substituting  (A4)  to  (A8)  into  (A3)  produces  the  prefactor  (i  for  a  fluid-elastic  interface. 

A. 2  The  Di rich  let  boundary  conditions 

The  precise  form  of  the  prefactor  J3  for  an  acoustic  field  incident  on  Dirichlet,  Neumann 
and  two-fluid  boundary  conditions,  and  for  an  electromagnetic  field  incident  on  a  perfect 
conductor  can  be  read  off  from  the  results  given  in  Appendix  C.l  of  reference  [6].  In 
particular,  the  result  corresponding  to  the  Dirichlet  boundary  conditions  is 

P  =  4 kizqiz  =  4 k0  sin  0m  sin  £out .  (A9) 
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